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ABSTRACT 

We present a study of the spatial distribution and kinematics of star-forming galaxies in 30 massive clus¬ 
ters at 0.15<z<0.30, combining wide-field Spitzer 24pm and GALEX NUV imaging with highly-complete 
spectroscopy of cluster members. The fraction ifsp) of star-forming cluster galaxies rises steadily with cluster¬ 
centric radius, increasing fivefold by 2 r 2 QQ, but remains well below field values even at 3 r 2 oo- This suppression 
of star formation at large radii cannot be reproduced by models in which star formation is quenched in in¬ 
falling field galaxies only once they pass within raoo of the cluster, but is consistent with some of them being 
first pre-processed within galaxy groups. Despite the increasing / 5 / 7 -radius trend, the surface density of star¬ 
forming galaxies actually declines steadily with radius, falling ^15x from the core to 2 r 2 oo- This requires 
star-formation to survive within recently accreted spirals for 2-3 Gyr to build up the apparent over-density of 
star-forming galaxies within clusters. The velocity dispersion profile of the star-forming galaxy population 
shows a sharp peak of 1.44 cr^ at 0 . 3 r 5 oo, and is 10-35% higher than that of the inactive cluster members at all 
cluster-centric radii, while their velocity distribution shows a flat, top-hat profile within rsoo. All of these results 
are consistent with star-forming cluster galaxies being an infalling population, but one that must also survive 
~0.5-2 Gyr beyond passing within r 2 oo- By comparing the observed distribution of star-forming galaxies in the 
stacked caustic diagram with predictions from the Millennium simulation, we obtain a best-fit model in which 
SERs decline exponentially on quenching time-scales of 1.73±0.25 Gyr upon accretion into the cluster. 

Keywords: galaxies; active — galaxies; clusters; general — galaxies; evolution — galaxies; stellar content 


1. INTRODUCTION 

The ability of galaxies to continuously form stars depends 
strongly on their global environment, with isolated central 
galaxies primarily evolving as star-forming spirals, while “red 
and dead” early-type galaxies completely dominate the cores 
of rich cl usters, producing the SE-density or SE-radius rela¬ 
tions (e.g. lKennicuttlfToOlLewis et alJl2002h . Various physi¬ 
cal mechanisms have been proposed over the years to remove 
(or consume) gas and quench star formation in spiral galaxies 
within massive clusters, such as ram-pressure or viscous strip¬ 
ping, st arvation, harassment or t i dal interactions (fo r reviews 
see e.g. iBoselli & Gavazzill2006tlHaines et al.ll2007h . 

Clusters and their member galaxies however do not ex¬ 
ist and evolve in isolation from the rest of the Universe. 
In ACDM models structure formation occurs hierarchi- 
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cally, meaning that as the most massive collapsed structures 
in the Universe, galaxy clusters for m latest and are also 
the most dynamically immature (e.g. iBovlan-Kolchin et alJ 
120091 iGao et a^]l2012^ . Preferentially residing at the nodes 
of the complex filamentary web, they continually accrete 
dark matter (DM) halos hosting individual ^L* galaxies 
(Mdm"^10'^Mq) or galaxy groups (Mdm~ 1 O'^“'"^M 0 ). The 
most massive clusters have on aver age doubled in mass since 
z~0.5 (iBovlan-Kolchin et al.ll200^ . while half of galaxies in 
local clusters have been accreted since z^O.4 (iBerrier et alJ 

HqoI. 


To correctly interpret the observed evolutionary and radial 
trends in cluster galaxy properties, it is thus fundamental to 
place them in this cosmological context whereby star-forming 
galaxies are being continually accreted into the clusters and 
transformed. Moreover, it is also vital to consider projection 
effects as many spectroscopic cluster members are actually 
infalling galaxies physically located outside the virial radius, 
and this contribution varies strongly with projected cluster¬ 
centric radius (rproj) and line-of-sight (LOS) velocity relative 
to the cluster redshift (v/o. 5 “(p))- This requires using cosmo¬ 
logical simulations containing one or more massive clusters, 
and following the orbits and merger histories of the galax¬ 
ies o r sub-halos which ar e accreted into the cluster over time 
(e.g. iMamon e t al.ll201 (ih. This approach ga i ned early sup- 

J iort when Balo gh. Navarro & Morrisl (120001) . iDiaferio et alJ 
200 ll) and Ellingson et al.l (120011) were able to reproduce the 
observed radial population gradients of star-forming galaxies 
in clusters in terms of galaxies on their first infall into the 
cluster. 

The caustic diagram, which plots vios-{v) versus rp^j, 
has been used to constrain the kinematics and ac¬ 
cretion epochs of different cluster galaxy populations. 
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as well as the constrain the masses, density profiles 


Moss & DickensI 119771: iBineseli et al.l 119871: iBiviano et al. 

1997 

,12013; Biviano & Kateerti 2004; Gill. Knebe & Gibson 

2005 

;lMahaia 

n et alJl201 It Hernandez-Eernandez et al.ll2014t 

Muzzin et all 

2014ll.laffeet alJl20L5l). 


This progress has permitted recent attempts to constrain 
the time-scales required to halt star formation in recently ac¬ 
creted cluster spirals, with results supporting gentle physi¬ 
cal mechanisms (e.g. starvation) that slo wly quench star- 


formation over a period of several Gvrs (iWol 

Fet all 120091 

von der Linden et alJ2010tlDe Lucia et al.ll2012t 

Wetzel et all 

20131). rather than more violent processes fe.e. mere- 


ers) that rapidly terminate star formation (although see 
e .g. iBalogh et al]r2004t iMcGee et al.ll201 it IWiiesinghe et alJ 


12012 ). 


In 


Haines et al.l (120091) we estimated the composite radial 


population gradients in the fraction of star-forming galaxies 
ifsp) in 22 massive clusters at 0.15<z<0.30 from the Local 
Cluster Substructure Survey (LoCuSSTHl based on panoramic 
Spitzer/MIPS 24pm data. A steady systematic increase in fsp 
with cluster-centric radius was observed out to ^r^nn, similar 
to those found previously ([Ellingson et al.ll2()()ll iLewis et aP 
I2002I IWeinmarin et al.ll2006l) . Bv comparison to galaxies in¬ 
falling and orbiting around massi ve clusters (M^nn^.lO ^^Mfr^) 
from the Millennium Simulation (ISpringel et al.ll20(CT . it was 
possible to approximately reproduce the radial population 
trends in the context of a simple infall model, in which star¬ 
forming field galaxies are accreted into the cluster and their 
star-formation rapidly quenched upon their first pericenter. 
The key limitation of this work was the lack of redshifts to 
identify cluster galaxies, such that we had to statistically ac¬ 
count for the contamination for field galaxy interlopers when 
estimating the fspif)- 

We have since completed ACReS (Arizona Cluster Red- 
shift Surve\Pl) which provides highly-complete spectroscopy 
of cluster members for all 30 clusters from Lo CuS!^ with 
wide-held Spitzer/MIPS data. With this data, iHaines et al.l 
(120131) found the specihc-SFRs of massive (Ad>10'°MQ) 
star-forming cluster galaxies within r 2 oo to be systematically 
28% lower than their counterparts in the held at hxed stellar 
mass and redshift, a difference signihcant at the 8.7tT level. 
This is the unambiguous signature of star formation in most 
(and possibly all) massive star-forming galaxies being slowly 
quenched upon accretion into clusters, and was best ht by 
models in which their star formation rates decline exponen¬ 
tially on quenching time-scales in the range 0.7-2.0Gyr. 

In this article we analyse the spatial distribution and kine¬ 
matics of star-forming galaxies within the same set of 30 
clusters, and by comparing with predictions from cosmolog¬ 
ical simulations, draw further independent constraints on the 
quenching time-scale. In particular, we determine the radial 
surface density prohle, S(r), of star-forming cluster galaxies 
and show that it declines steadily with radius, falling ^15 x 
from the core to 2r2oo- We show that this simple obser¬ 
vation provides powerful constraints for how long massive 
star-forming galaxies are able to continue forming stars once 
they are accreted into rich clusters, quickly ruling out models 
in which star-formation is rapidly halted in infalling spirals 
when they pass within r 2 oo- We also re-examine the radial 
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population gradients of star-forming galaxies (/^/r-radius re¬ 
lation) out to 3 r 2 oo, where we find a shortfall of star-forming 
galaxies in comparison to the coeval field population that can¬ 
not be easily explained by purely cluster-related quenching 
mechanisms, indicating a need for galaxies being first pre- 
processed within infalling galaxy groups. 

In ^we present our observational data, and in ^the main 
results. In ®we follow the infall and orbits of galaxies in the 
vicinity of massive galaxy clusters from the Millennium simu¬ 
lation, to predict their spatial distributions and kinematics as a 
function of accretion epoch. These model predictions are then 
compared to observations in |5] We discuss the resultant con¬ 
straints on the time-scales required to quench star formation 
in recently accreted galaxies and the need for pre-processing 
in ^and summarize in ^ Throughout we assume Hm= 0.3, 
r2A=0.7 and Ho=70kms“'Mpc“'. 

2. DATA 

LoCuSS is a multi-wavelength sur vey of X-ray lumi- 
nous galaxy clusters at 0.15<z<0.3 (ISmith et al.l 12010a^ 
drawn from the ROSA T All Sky Survey cluster catalogs 
(iBohringer et al.l 120041) . The first 30 clusters from our 
survey benefit from a particularly rich dataset, including : 
Subaru/Suprime-Cam optical imaging (lOkabe et alJ l2010i) . 
Spitzer/MIPS 24pm maps, //ersc/ie//PACS-)-SPIRE 100- 
500/rm maps, Chandra and/or XMM X-ray data, GALEX UV 
data, and near-infrared (NIR) imaging. All of these data em¬ 
brace at least 25'x25' fields-of-view centered on each clus¬ 
ter, and thus probe them out to 1-2 virial ra dii (IHaines et alJ 
l2010t|Pereira et af]l2010t [Smith et ani2010bh . These 30 clus¬ 
ters were selected from the parent sample simply on the basis 
of being observabl e by Subaru on the nights allocated to us 
(lOkabe et al.ll2010h . and should therefore not suffer any gross 
biases towards (f or example) cool co re clusters, merging clus¬ 
ters etc. Indeed, lOkabe et alJ (120101) show that the sample is 
statistically indistinguishable from a volume-limited sample. 

2.1. Chandra/XMMX-ray imaging 

All but two (Abell 291, Abell 2345) of the 30 clusters 
have available deep Chandra data (f„p=9-120ksec). De- 
projected dark matter densities, gas densities and gas tem¬ 
perature profiles for each cluster were derived by fitting 
the ph enomenological cluster models of lAscasibar & Diegol 
( 20 ()8h to a series of annular spectra extracted for each cluster 
(ISanderson & PonmarJl20lH) . The best-fitting cluster models 
were then used to estimate r^oo, the radius enclosing a mean 
overdensity of 500 with respec t to the critical density of the 
Universe at the cluster redshift (ISanderson et al.ll2009[ f _ 

The rsoo value for Abell 689 is taken from iGiles et alJ 
dloTl) . as they separated the extended cluster X-ray emission 
from the central BE Lac, and for the clusters lacking Chan- 
dra data, the rsnn va lues are taken from the XMM analysis of 
iMartino et alJ Eor the 23 clusters in common with the 

joint Cha ndra-XMM analysis o f the LoCuSS high-Lx cluster 
sample of IMartino et al.l (120141) . there is good consistency of 
the cluster radii with (r5oo,/ra,«ei/i-5oo,Mart/«o)=E003±0.034. 

The Chandr a data were a l so use d to identify X-ray AGN 
as described in IHaines et al.l (l2012h . The survey limit of six 
broad (0.3-7keV) X-ray counts results in on-axis sensitivity 
limits of Lx<L0x 10"^^ergs“* for X-ra y AGN at the clust er 
redshift for all 28 systems (Table 1 from iHaines et al.ll2012l) . 

Deep XMM data was available for 23 systems, allow¬ 
ing other groups and clusters in the region to be identi¬ 
fied. Each 0.5-2 keV image is decomposed into unresolved 
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Figure 1. Radial coverage of the Spitz^r/MlPS 24fim data. Each curve shows 
the fraction of cluster members covered by our 24/im images (25' x 25' field- 
of-view) as a function of cluster-centric radius in units of r^QQ, for clusters 
in the following redshift ranges: 0.15<z<0.20 {long-dashed blue curx’es)', 
0.20<z<0.25 {dot-dashed green curves)', 0.25<z<0.30 {short-dashed red 
curves)', 0.15<z<0.30 {solid black curx’es). 

and extended em i ssion, following the wavelet technique of 
iFinoguenov et alJ (120091) . For each extended source, we at¬ 
tempt to identify the redshift of its associated group/cluster 
by examining the Subaru optical images for likely BCGs near 
the center of the X-ray emission and/or groups of galaxies 
with similar redshifts from ACReS within the X-ray contours. 

2.2. Mid-infrared Observations 

All 30 clusters were observed at 24^m with MIPS 
(iRieke et al.l 12004) on board the Spitzer Space Telescop^ 
(PID: 40872; PI: G.P Smith). T he resulting 24rtm mo- 
saics were analysed with SExtractor (iBertin & Arnoutsll996h 
producing catalogs which are on average 90% complete to 
400/rJy. 

Each cluster was observed across a fixed 25'x25' field-of- 
view, resulting in the clusters being covered out to different 
cluster-centric radii in units of r^oo, depending on their red¬ 
shift and rsoo radius, as well as the orientation of the Spitzer 
images. Eigure [1] shows that we probe out to larger cluster¬ 
centric radii for the highest redshift clusters (0.25<z<0.30) 
than those in our lowest redshift bin (0.15<z<0.20). Averag¬ 
ing over the full redshift range (0.15-0.30; solid black curve), 
our 24/im coverage is essentially complete out to r 2 oo, falling 
to ^45% at 2 r 2 oo, based on the conversion rsno = 0.66 r 2 oo 
(iSanderson & Ponmanll2003h . 

2.3. UV, optical and near-infrared data 

Wide-held J- and X-band near-infrared imaging was ob¬ 
tained for all 30 clusters using either WECAM on the 3.8- 
m United Kingdom Infrared Telescope (UKIRT0 (52'x52' 
held-of-view; 26/30 clusters) or NEWFIRM on the 4.0-m 
Mayall telescope at Kitt PealQ (27'x27'; 4/30 clusters), in 
each case reaching depths of X^19, y~21. 

Wide-held deep UV imaging from the Galaxy Evolution 
Explorer (GALEX) satellite was obtained for 26/30 clusters, 

Tliis worlt is based in part on observations made witli the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA (contract 1407). 

UKIRT is operated by the Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom. 

Kitt Peak National Observatory, National Optical Astronomy Observa¬ 
tory, which is operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the National Science 
Foundation. 


primarily through the Cycle 4 (GI4-090; PI. G.P. Smith) and 
Cycle 6 (GI6-046; PI. S. Moran) Guest Investigator Pro¬ 
grams. The Cycle 4 program provided far-ultraviolet (FUV) 
and near-ultraviolet (NUV) imaging for 14 clusters (f„p=3.2- 
13.6ksec), while comparable FUVh-NUV data were obtained 
for 7 more clusters (f„p=3.4-29.0ksec) from the GALEX sci¬ 
ence archive. Sixteen clusters were observed in Cycle 6, 
including 8 systems not previously observed, but this pro¬ 
vided only deep NUV imaging (f(,;ip=2.9-36.3ksec), the oper¬ 
ations of the FUV camera having previously been suspended. 
The GALEX instrument has a circular field of view of radius 
0.55 deg, ensuring full ultraviolet coverage for galaxies in our 
near-infrared WECAM fields. The total NUV exposure times 
and 5tT magnitude lim i ts, afte r correcting for Galactic extinc¬ 
tion as in IWvder et al.l (120071) . are shown in Table[T] 

Optical photometry in the ugriz bands were taken from 
the Sloan Digital Sky Survey, using the dereddened SDSS 
model magnitudes. Twenty-six of the 30 clusters lie within 
the DR-10 footprint, while 23 have both SDSS ugriz and deep 
GALEX NUV photometry, allowing star-forming galaxies to 
be identified from their blue NUV-r colors. 

2.4. MMT/Hectospec spectroscopy 

We have recently completed ACReS (the Arizona Cluster 
Redshift Survey; Pereira et al. 2015 in preparation) a long¬ 
term spectroscopic programme to observe our sample of 30 
galaxy clusters with MMT/Hectospec. Target galaxies are 
primarily X-band selected down to a limit of m|.(zc;)+1.5 or 
fainter (depending on the number of targets produced), to 
produce an approximately stellar mass-limited sample down 
to XI^2xlO'®M0. Higher priorities are given to target 
galaxies also detected at 24/im to obtain a virtually com¬ 
plete census of obscured star formation in the cluster popula¬ 
tion. Further details of the survey aims and targetting strategy 
are given in iHaines et al.l (120131) . Eleven of our 30 clusters 
were also observe d by the Hectospec Cluster Survey (HeCS; 
iRines et dl 120131) . providing redshifts for an additional 971 
cluster members. Redshifts for a further 112, 92 and 49 mem¬ 
bers of clusters 1 ^11720. l-b2638. Abell 383 and Abell 209 
are included fro mlOwers et akl (1201 Ih . lGeller et ai1(^2014^ and 
iMercurio et al.l (120031) respectively. Table [T| lists the number 
of spectroscopic members for each cluster, with the contribu¬ 
tions taken from other published surveys indicated in paren¬ 
theses, giving us a grand total of 10 950 cluster members with 
redshifts. Averaging over all 30 systems, we achieve spectro¬ 
scopic completeness levels of 66% foxM k<— 21>.\0 {M’f-vl .5) 
cluster galaxies across the full WFCAM/NEWFIRM fields, 
rising to 80% for those galaxies with Spitzer coverage and 
96% for those detected at 24/rm. 

The likelihood that a given galaxy was targetted for spec¬ 
troscopy depends strongly on both its location with re¬ 
spect to the cluster center as well as its photometric prop¬ 
erties (X -band magnit u de, J- K color, 24/im flux), as de¬ 
tailed in IHaines et all (120131) . To account for this, each 
galaxy is weighted by the inverse probability of it having 
being observed spec troscopically, following the approach of 
iNorberg et al.1 (l2()02h . 

2.5. Identification of cluster members and field galaxy 
samples 

Members of each cluster are identified from the redshift 
versus projected cluster-centric radius plot as lying within 
the “trumpef’-shaped caustic profile expected for galaxies in¬ 
falling and subsequently orbiting within a massive virialized 
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Cluster 

(z> 

Nz (lit) 

^•500 

Mjoo 


texp(NUV) 

mAB(NUV) 

Name 



(Mpc) 

(10*''Mq) 

(km ) 

(sec) 

(S/N=5) 

Abell 68 

0.2510 

194 (0) 

0.955 

3.193 

1186f»g 

6731 

22.90 

Abell 115 

0.1919 

213 (36) 

1.304 

7.628 

1219)^2 

16217 

23.36 

Abell 209 

0.2092 

393 (49)^ 

1.230 

6.519 

1369)“ 

11616 

23.16 

Abell 267 

0.2289 

230 (139)“ 

0.994 

3.515 

1045+^2 

11268 

23.32 

Abell 291 

0.1955 

126 (0) 

0.868" 

2.259 

704)*° 

6059 

23.16 

Abell 383 

0.1887 

266 (92)'' 

1.049 

3.958 

950)“ 

6117 

23.16 

Abell 586 

0.1707 

247 (21) 

1.150 

5.117 

933)“ 

— 

— 

Abell 611 

0.2864 

297 (7) 

1.372 

9.847 

1039)“ 

9928 

23.01 

Abell 665 

0.1827 

359 (31) 

1.381 

8.975 

1227+g 

9503 

23.26 

Abell 689 

0.2776 

338 (153)“ 

1.126^ 

5.390 

721)“ 

— 

— 

Abell 697 

0.2821 

486 (141)“ 

1.505 

12.93 

1268)“ 

36858 

23.89 

Abell 963 

0.2043 

466 (50)" 

1.275 

7.226 

1119)^° 

29043 

23.96 

Abell 1689 

0.1851 

857 (416)“ 

1.501 

11.55 

1541)1° 

7716 

23.34 

Abell 1758 

0.2775 

471 (50)" 

1.376 

9.835 

1442)“ 

20977 

23.88 

Abell 1763 

0.2323 

423 (126)“ 

1.220 

6.522 

1358)“ 

13589 

23.80 

Abell 1835 

0.2524 

1083 (608)“ 

1.589 

14.73 

1485)“ 

21827 

23.50 

Abell 1914 

0.1671 

454 (65)" 

1.560 

12.73 

1055)11 

3804 

23.17 

Abell 2218 

0.1733 

342 (49) 

1.258 

6.716 

1245)12 

14617 

23.73 

Abell 2219 

0.2257 

628 (297)“ 

1.494 

11.90 

1332)“ 

9886 

23.42 

Abell 2345 

0.1781 

405 (39) 

1.249" 

6.607 

1000)11 

9864 

23.49 

Abell 2390 

0.2291 

517 (140) 

1.503 

12.16 

1372)“ 

6819 

22.35 

Abell 2485 

0.2476 

196 (0) 

0.830 

2.088 

799)54 

— 

— 

RXJ0142.0-I-2131 

0.2771 

204 (15) 

1.136 

5.531 

1123)“ 

— 

— 

RXJ1720.1+2638 

0.1599 

473 (114)* 

1.530 

11.92 

938)“ 

3074 

22.67 

RXJ2129.6+0005 

0.2337 

334 (78)" 

1.227 

6.648 

879)“ 

36323 

23.72 

ZwClO 104.4+0048 (Z348) 

0.2526 

185 (1) 

0.760" 

1.613 

806+11 

14435 

23.46 

ZwC10823.2+0425 (Z1693) 

0.2261 

337 (4) 

1.050 

4.130 

671^0 

14314 

23.41 

ZwC10839.9+2937 (Z1883) 

0.1931 

173 (3) 

1.107 

4.674 

834)“ 

21250 

23.39 

ZwC10857.9+2107 (Z2089) 

0.2344 

147 (0) 

1.024 

3.866 

815)“ 

20225 

23.68 

ZwC11454.8+2233 (Z7160) 

0.2565 

157(1) 

1.128 

5.294 

988)“ 

19850 

23.58 


Table 1 

The cluster sample. Column (1) Cluster name; col. (2) Mean redshift of cluster members; col. (3) Total number of spectroscopic cluster members (contribution 


taken from the literature, including new redshifts from IRines et ^ <20131) ^. lOwers et al.l d2011h ^. lMercurio et al]d2003l) ^ and IGeller et al.1 col. (4) 

radius rsoo in Mpc. ^From lMartino et ^d2014l) . ^From Giles et al. (2012); col. (5) Cluster mass M 500 in IO^'^Mq ; col. (6) Velocity dispersion of cluster 
members within r 2 oo', col. (7) Total exposure time of GALEX NUV images; col. (8) NUV magnitude limit for SNR = 5 


structure (iDiinner et alJ I2007h . For most systems, there is 
a strong contrast i n phase-space density fr om inside to out¬ 
side these caustics (iRines & Diaferioll2006l) . making their vi¬ 
sual identification relatively simple. The central redshift (z) 
and velocity dispersion of each cluster (Table [Til are iter- 
ative ly measured for member galaxies within r 2 oo (estimated 
as iniFinn et al.l2 008h. using the biweight scale estimator (Sbi', 
iBeers, Flvnn & Gebhardlill99(lh . with uncertainties estimated 
using bootstrap resampling. 

The field galaxy sample was taken from the same dataset 
as the primary cluster galaxy sample, but were located in nar¬ 
row redshift ranges on either side of the cluster, for which 
our spectroscopic survey ACR eS should still be c omplete to 
M|--I- 1.5 (for full details see iHaines et al.ll2013h . Overall, 
1398 coeval (0.15<z<0.30) field galaxies withMA:<M^-l-1.5 
and 24pm coverage were identified within these narrow red¬ 
shift slices either side of the clusters (699 in front, 699 
behind), after excluding regions where other X-ray galaxy 
groups had been previously detected from our XMM data. 

2.6. Stellar masses and SFRs 

Rest-frame UV-optical colors and ab solute magnitudes 
were d etermined using the k-corrections of IChilingarian et Tl 
(l201(]h . Stellar masses (M) were estimated from the X-band 
luminosities using the linear relation between X-b and stellar 
mass-t o-light ratio and rest-frame g - i color from iBell et al.l 
(120031) . adjusted by -0.15 dex to be valid for a lKrounal (l2002h 


IMF. Where SDSS photometry was unavailable we classified 
the galaxy as being either star-forming or passive according 
to whether it was 24pm or NUV detected or not, and adopted 
appropriate M/L ratios. 

For each 24pm-detected galaxy with known redshift, its in¬ 
trinsic bolometric luminosity (Lr/s) and rest-frame 24pm lu¬ 
minosity is estimated by comparison of its 24pm flux to the 
luminosit y-dependent t empla te infrared spectral energy distri¬ 
butions of iRieke et al.l (l2009h . The latter is then converted t o 
an obscured SFR using the calibration of iRieke et al.l (l2009h 

SFR/«(MQyr“') = 7.8x 10-“’L(24pm,LQ) (1) 

which is valid for either a iKroupal (l2002h or IChabrieil (120031) 
IMF. Our Spitzer data should be sensitive to galaxies with 
ongoing obscured star formation occuring at rates down to 
2 .OM 0 yr“' in our most distant clusters (z~0.28). 

Local quiescent early-type galaxies are known to emit in 
the mid-infrared at levels m uch higher than expe cted from 
photospheric emission alone (IClemens et al.ll200^ . This ex¬ 
cess at 10-30pm has been shown to be due to silicate emis¬ 
sion from the dusty circumstellar envel opes of mass-losing 
evolved AGB stars (iBressan et al.ll2006l) rather than residual 
ongoing star formation. We may thus worry that some of our 
cluster galaxies may be mistakenly classed as star-forming 
due to 24pm emission coming from TP-AGB stars. 

The spectral energy distributions (SEDs) of evolved stel- 
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Figure 2. Radial population gradients for MIR-selected star forming galax¬ 
ies from our stacked sample of 30 clusters. Red symbols show the fraction of 
massive {Ai > 2 .Ox IO^^Mq) cluster galaxies with obscured star-fonnation 
at rates SFR/^ > 2.0Mq yr ^ as a function of projected cluster-centric radius 
(fproy/fsoo)- The error bars indicate the un certainties deriv ed from binomial 
statistics calculated using the formulae of lGehreli 419861) . Each radial bin 
contains 400 cluster galaxies. The blue horizontal line indicates the corre¬ 
sponding fraction of field galaxies (A4>2.0x 10*® Mq; 0.15<z<0.30) with 
SFR/^i > 2 . 0 Mq yr“* and its Itr confidence limits (shaded region). 

lar populations including emission from dusty circumstel- 
lar envelopes peak at 10 - 20 rtm, but th en drop rapidly 
at longer wavelengths (iPiovan et al.l l2003h . and so galax¬ 
ies whose 24/rm emission is due to TP-AGB stars should 
not be detected in our Herschel/PACS data, unlike nor¬ 
mal star-form ing galaxies who se infrared SEDs peak at 70- 
170pm (e.g. iDale et all 120121) . For the 11 nearest clus¬ 
ters (0.15<z<0.20) in our sample, >98% of galaxies with 
SFR/s> 2 .OM 0 yr“* and HerschelfPACS coverage were also 
detected at 100pm, while >99% show clear Ha emission 
in our ACReS MMT/Hectospec spectra, indicating that their 
24pm emission is indeed due to ongoing star formation. 

In a comparable Spitzer/MIPS analysis of 814 galaxies 
in the Shapley supercluster at z=0.048, sensitive to much 
lower obscu r ed star formation rates (SFR/R^O.O 5 M 0 yr“*) 
iHaines et akl (1201IH) did find a significant population of qui¬ 
escent (based upon a lack of Ha emission) cluster galaxies de¬ 
tected at 24pm, but none with 24pm luminosities that would 
convert to an obscured SFR above 0.5 M 0 yr“*. They also ob¬ 
tained a tight correlation (0.22 dex) between the 24pm and the 
1.4 GHz radio luminosities for star-forming cluster galaxies, 
down to SFR/R~l.OM 0 yr“*, consistent with both the mid- 
infrared and radio emission being due to star formation. 

3. RESULTS 

3.1. Radial population gradients 

Figure|2]shows the fraction of massive (yVl>2.0x 1 O*°M 0 ) 
cluster galaxies with obscured star formation occurring at 
rates SFR/s > 2 .OM 0 yr“* (fsp) as a function of projected 
cluster-centric radius in units of rsoo, across the full sample 
of 30 clusters. The fractions only consider cluster members 
covered by the Spitzer 24pm maps, and exclude X-ray AGN 
and QSOs, as their 24pm emission is usually dominated by 



Figure 3. X’-band luminosity functions (LFs) for cluster galaxies within rsoo 
(red solid curx’e), galaxies from the infall regions (rproj>ri 0 Q\ green solid 
curve) and the coeval field population (light-blue solid curx’e). The latter two 
are normalized to contain the same total number of <—23.1 galaxies as the 
cluster galaxy population, and the points slightly shifted horizontally, to per¬ 
mit easie r compar i son of the LFs. Enor bars indicate Poissonian uncertaintes 
based on IGehreT3 419861) . The corresponding A'-band LFs for star-forming 
galaxies with SFR;^> 2 .OM 0 yr“* are indicated by dashed curves. 

dust heated by the active nucleus rather than star formation 
(e.g. Xu et al. 20 15). BCGs are als o excluded due to their 
unique evolutions (iFin & Mohiil2004l) . and the direct link be¬ 
tween B CG activity and the presence of coolin g flows within 
clusters (iSmith et alJl20103:lRawle et al.ll2012h . 

The fraction of obscured star-forming galaxies increases 
steadily with cluster-centric radius from fgf'^O.OA in the clus¬ 
ter core to fsf^Q.23 at 3 . 0 r 5 oo ( 1 . 9 r 2 oo)- However even at 
these large radii the fgf remains well (^1/3) below that seen 
in coeval field galaxies (/y^ =0.33±0.01; blue line). A simple 
linear extrapolation of the observed trend suggests that the fgp 
should reach that of the field galaxy population at ^ 4 . 5 r 5 oo- 
However the limited extents of our Spitzer 24pm maps mean 
that we cannot establish whether this occurs or not. 

The fraction of star-forming cluster galaxies evolves very 
rapid ly at these redshifts, with/y^(x(l-l-z)^ ®^* ' (IHaines et alJ 
lMT3h . It is thus vital to ensure that this shortfall in star¬ 
forming cluster galaxies at large radii with respect to field 
values is not produced by a redshift bias between the two sam¬ 
ples. This is certainly not the case here, as the mean redshifts 
of each radial bin for the cluster populations all lie in the range 
0.217-0.241, while for the field population (z) = 0.225. There 
is a marginal redshift bias within our cluster sample, as the 
outer two radial bins have (z) = 0.241, while the remaining 
bins all have mean redshifts in the range 0.217-0.231. This 
is due to our Spitzer data providing wider radial coverage (in 
terms of ryoo) for the higher redshift systems, but it is likely 
only a marginal effect, artificially increasing the outer two fsp 
by < 10 % (or ^ 0.02 in the figure). 

A second possible explanation for the lower fgf among 
cluster galaxies would be if they were more massive on aver¬ 
age than the field galaxy comparison sample, as fgf is known 
to d ecline with increasi ng stellar mass at fixed galaxy density 
(e.g. IHaines et al.ir2007t) . However, the mean stellar masses of 
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Figure 4. Comparison of radial population gradients for star-forming clus¬ 
ter galaxies selected in the ultraviolet and mid-infrared. Magenta squares 
show the fraction of massive (A1>2.0x IO^^Mq) cluster galaxies with blue 
rest-frame UV-optical colors, {NUV-r)Q_o<4.5, indicating ongoing unob¬ 
scured star-formation. Each radial bin contains 300 cluster galaxies. The 
magenta dashed line indicates the corresponding fraction of field galaxies 
(A^>2.0xl0i”Mo; 0.15<z<0.30) with (A^f/y-r)o.o<4.5 and its la confi¬ 
dence levels (magenta shaded region). Grey symbols show the corresponding 
fractions of cluster galaxies with SFR//? > 2.0MQyr“^ taken from Fig. [2| 
scaled via the right-hand axis, such that the fraction of field galaxies with 
SFRjr > 2.0MQyr * (blue dashed line) coincides with the corresponding 
fraction selected in the NUV. 


cluster galaxies for each radial bin in Fig.l^are always within 
0.04 dex of that of the coeval field population. Moreover, the 
/T-band luminosity functions of cluster galaxies within 1-500 
(excluding BCGs), those in the infall regions (rpioj>r 50 Q), and 
the field galaxy samples are all indistinguishable (Fig. (3). We 
also note that due to the fixed SFR lower limit of 2.0 Mq yr“' 
used to define a star-forming galaxy, the fsf do not vary much 
with stellar mass (A'-band luminosity). We therefore exclude 
secular quenching due to increased stellar masses among the 
cluster galaxy population as being responsible for their lower 
fsf (at all radii in Fig. ^ with respect to the field. The pro¬ 
gressive suppression in star-formation in moving from field 
to infall regions, and on to cluster environments is seen at all 
stellar masses (compare dashed curves in Fig. [3. 

The GALEX NUV data provides a complementary means of 
identifying star forming galaxies from their ultraviolet emis¬ 
sion, and an opportunity to measure the SF-radius relation out 
to larger cluster-centric radii. Figure 0] shows the composite 
radial population gradient {magenta squares) in the fraction 
of cluster galaxies having blue rest-frame UV-optical col¬ 
ors {NUV—r)o,o<4.5, for the 23 clusters with deep GALEX 
NUV data and SDSS ugriz photometry. The (Nt/V-r)o.o=4.5 
color limit lies in th e middle of the UV-optical “green valley” 
(IWvder et al.ll200'^ . and allows passively-evolving galaxies 
to be efficiently excluded w ithout losing dusty star-forming 
galaxies due to reddening (iHaines et al.l l200a) . Optically- 
quiescent early-type galaxies with residual (or “rejuvenated”) 
extended star formation in th e form of rings or s piral arms 
should also still be recovered (ISalim & RichlUOllih . 

As before, the fraction of star-forming cluster galaxies in¬ 
creases steadily with cluster-centric radius from11 in 



clustercentric radius (rproj/rsoo) 

Figures. Composite galaxy surface density profile for all 
galaxies (BCGs excluded) from our sample of 30 0.15<z<0.30 clusters, as 
a function of projected cluster-centric radius (rproy/rsoo)- Each symbol con¬ 
sists of 200 spectroscopic cluster members. The blue symbols show the corre¬ 
sponding observed surface density S(r) of Ma:<M^+ 1.5 galaxies with active 
obscured star-formation at rates SFR//?>2.0MQyr“^ from the same ensem¬ 
ble cluster. The magenta squares show the radial profile for unobscured star¬ 
forming galaxies with (A^f/y-r)o.o<4.5 from the 22 clusters that have SDSS 
ugriz and deep GALEX NUV photometry. The en'or bars assume Poisson 
statistics. The corresponding best-fit NFW profiles are shown by the red, 
blue and magenta curves. The dot-dashed black curve shows the predicted 
surface density profiles for all “spectroscopic” member galaxies of the 75 
most massive clusters from the Millennium simulation, scaled to fit our ob¬ 
served ensemble surface density profile of all galaxies. 

the cluster core to /sf~0.44 at 4-5 1-500 (~3 r 2 oo)- The fraction 
of star-forming cluster galaxies remains significantly below 
that seen in coeval field galaxies (/5/r=0.53±0.01; magenta 
dashed line), even out at ~5r5oo, the trend appearing to flat¬ 
ten off rather than continue upwards to field values. Again the 
two samples are confirmed to be coeval: each radial bin for 
the cluster population has a mean redshift in the range 0.217- 
0.236, while for field galaxies (z) = 0.226. 

To allow comparison between the NUV-based and 24/im- 
based SF-radius relations, the latter {gray points) is replotted 
from Fig. |2] after adjusting its vertical scale {right hand axis) 
to ensure that fractions of IR-selected and UV-selected star¬ 
forming field galaxies coincide on the plot. While both re¬ 
lations show the same steadily increasing trends with radius, 
the NUV-based SF-radius relation consistently lies above the 
re-scaled IR-based relation. 

Looking at the IR-based SF-radius relation, it is tempting 
to suggest that the fraction of star-forming galaxies even falls 
to zero at the cluster core. However, we confirm that this is 
not the case. Further splitting each radial bin into four, the 
fsF,iR never fall below a floor value of ~3-5% in the cluster 
core. Similarly, the NUV-based SF-radius relation never falls 
below ~8-15% in the cluster core, when the radial bins are 
further sub-divided. A residual population of cluster galaxies 
with ongoing star formation exists at all radii. 

3.2. Radial surface density profiles 

The spatial distribution of galaxies within clusters provide 
key constraints on the primary epoch at which they were ac¬ 
creted into the system, as well as the effects of continual clus¬ 
ter mass growth on their orbital parameters and tidal and ram- 
pressure stripping on their stellar masses. Figure|3shows the 
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Figure 6. Stacked observed phase-space diagram, {^ios~{'^))l versus 
^pro ;/^5005 of galaxies combining all 30 clusters in our sample. Each gray 
solid dot represents a spectroscopic cluster member, while open points indi¬ 
cate non-cluster galaxies. Only those galaxies covered by our Spitzer data ai'e 
plotted. Star-forming galaxies detected at 24/im are indicated by larger sym¬ 
bols, colored according to their specific-SFRs from red (sSFR<10“^^ yuM 
to blue {sSFR>10“^^ yr“^). The blue curve and error bars indicates the l-cr 
velocity dispersion profile of the 24/im-detected cluster members. The un¬ 
certainty in each a{r) value is estimated by bootstrap resampling the galaxies 
in that radial bin. The black curve shows the corresponding radial profile for 
the remaining inactive cluster members. 

radial distribution of all Mk < M^+1.5 galaxies in the stacked 
cluster out to ^4r5oo, excluding BCGs (redpoints). 

Large numerical DM simulations have found that DM ha- 
los are well describ ed by a “universal” 2-parameter (NFW; 
iNavarro et al.lll997h density profile from scales of lOkpc out 
to lOMpc (lFrenketal.lfT9M iGao et al.llMTl) . The NFW 
profile is characterized by a scale radius rs=r 2 oo/c, where c 
is the concentration parameter. The three-dimensional den¬ 
sity profile is given by p(x) cx x~^(l+x)~^, where x=rlrs, 
and d\ogp/d\ogr = -2 at r=rs. The NFW model has been 
shown to provide excellen t fits to stacked tang ential shear pro¬ 
files of massive clusters (lOkabe et al.ll2013h and the distri¬ 
bution of cluster galaxies jLin, Mohr & Stanfordll2004ll . The 
surface density profile, S(r), of cluster galaxies in our en¬ 
semble cluster sample can be well described by a projected 
NFW profile with a concentration parameter Cg=3.01±0.16 
(re d curve), consistent with the C p=2.90±0.22 value obtained 
by iLin, M ohr & StanfordI (1200^ or Ci,=2.7zfc 0.7 obtained by 
iBudzvnski et al.l (1201 2h . while iMuzzin et akl (120071) obtained 
a higher concentration of Cg=4.13±0.57 for the stacked K- 
band number density profile of galaxies from 15 clusters at 
0.2 < z < 0.55. The Cg value of 3.01 obtained here is still 
significantly lower than the cn/L.200=4.22j^-3g value for the 
concentration of th e overall mass distribution obtained by 
lOkabe et al.l (1201 3h from their weak lensing analysis of 50 
z^O.2 clusters, including many of the systems in our sample. 

The surface density, E(r), of massive (Mk<M'^+1.5) clus¬ 
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Figure 7. Stacked relative LOS velocity distributions of inactive galaxies 
{red solid histograms), 24/im-detected {blue striped histogram) and UV- 
selected {NUV — r < 4.5; magenta striped histogram) star-forming galaxies 
from our ensemble cluster, in four bins of projected radius. The standard de¬ 
viations (in units of cFp,) and kurtosis values of each distribution are indicated. 

ter galaxies with obscured SFR/s>2.0MQyr“* (blue points) 
declines steadily with radius from ^7 galaxies r^^ cluster"' 
in the cluster cores to ^0.4 by ^3r5oo. These S(r) include 
corrections for spectroscopic incompleteness (§ 12.41 a nd the 
radial variation in coverage by the 24pm images (§ 12.21) . Even 
though the fraction of star-forming galaxies is falling to close 
to zero in cluster cores, clusters mark over-densities in the 
spatial distribution of star-forming galaxies in the plane of 
the sky. This is simply due to the cuspy density profile of 
the global cluster galaxy population more than compensat¬ 
ing for the steady decrease in fsr when approaching the clus¬ 
ter core. The surface density of Mk<M^+\.5 cluster galax¬ 
ies with (NUV - r)ox) < 4.5 (magenta squares) shows a sim¬ 
ilar radial profile, but is margina lly s teeper at small radii 
(rproj ^ 0.2r5oo)- We will show in § l5.2l that these steadily de¬ 
clining trends in S(r) imply that star formation must survive 
within recently accreted spirals for several Gyr to build up the 
apparent over-densities of star-forming cluster galaxies. 

3.3. Dynamical analysis of star-forming cluster galaxies 

Figure |6] shows the stacked caustic diagram of all 30 clus¬ 
ters, in which the projected radius of each cluster member is 
normalized by the Chandra-based rsoo of that cluster, and the 
LOS velocities are scaled in units of The 24pm-detected 
star-forming galaxies (colored points) do not have the same 
spatial distribution within the caustic diagram as the remain¬ 
ing inactive cluster galaxy population (gray points). They 
appear to preferentially lie along the caustics, indicative of 
an infalling population. They also show a concentration at 
-0-1.5 r 5 oo, covering the full velocity range within the 
caustics. The fall off in numbers towards larger radii is a se¬ 
lection effect due to the decline in Spitzer 24/im coverage be¬ 
yond 2r5oo. Star-forming galaxies do not entirely avoid the 
central region with low cluster-centric radii (rpro/<0.4r5oo) 
and relative L( 3S velocities (|Ati/ cr,J<,0.80), as X-ray AGN 
appear to do (iHaines et al.l 12012) . but their frequency cer¬ 
tainly drops off here, in marked contrast to the inactive cluster 
galaxy population. Some of the star-forming galaxies in these 
central regions of phase-space will likely appear here due to 
projection effects, being located along the line-of-sight of the 
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cluster but physically still well outside r 2 oo (see Fig.l22]). 

The velocity dispersion of the star-forming cluster galaxy 
population {blue curve) is 10-35% higher than that of the in¬ 
active cluster galaxies {gray curve) at all radii. Averaging over 
all galaxies within 1-500 ( 21 - 500 ), star-forming galaxies have ab¬ 
solute LOS velocity offsets which are 26.4% (24.9%) higher 
than their passive counterparts at the same cluster-centric ra¬ 
dius, a result significant at the S.Ocr (10.Ter) level. 

The same trends are obtained when selecting star-forming 
cluster galaxies according to their rest-frame NUV - r color. 
Their velocity dispersion remains 10-35% higher than that of 
inactive cluster galaxies out to 2.5-3 r 5 oo, although at larger 
radii they become indistinguishable. 

The relative LOS velocity distributions of star-forming 
(both 24/im-detected and UV-selected) and inactive cluster 
galaxies for the stacked LoCuSS cluster sample are shown in 
Figure|7] in four bins of projected cluster-centric radius. The 
LOS velocity distributions of inactive galaxies (not detected 
at 24/im and having {NUV-r)Q,o > 4.5) can be approximately 
described as a Gaussian at all radii. The distributions of star¬ 
forming galaxies in the two inner radial bins {rproj<r 5 Qo) in¬ 
stead appear more consistent with a flat, top-hat profile than a 
Gaussian, including an relative excess of star-forming galax¬ 
ies at |Aj//(Ti/| > 1.2 in comparison to the inactive popula¬ 
tion. At 0 . 5 -LOr 5 oo there is even marginal evidence a cen¬ 
tral dip in the LOS velocity distribution of star-forming galax¬ 
ies. Thekurtosis ,72 = 

of the LOS velocity distributions of star-forming galaxies 
( 72 ~“ 0 . 81 ) is significantly lower than that expected for a 
Gaussian distribution ( 72 = 0 . 0 ) at >5 ct level in both inner ra¬ 
dial bins, and closer to expectations for a uniform top-hat 
distribution ( 72 =- L2). At large cluster-centric radii (2.0- 
4.0 r 5 oo), the velocity distribution of star-forming galaxies is 
consistent with a Gaussian function, and is almost indistin¬ 
guishable from that of the inactive population, albeit with a 
marginally (^9%) higher dispersion. 

To further delineate the connection between the kinematics 
of cluster galaxies and their ability to form stars, the sample 
of 24/rm-detected cluster galaxies within 2 . 0 r 5 oo is split into 
four bins of specifle-SFR and their mean cluster-centric radii 
{fproj /rsoo), and LOS velocity dispersions cfios{<2rsQo) /(Jv, 
compared in Figured Kinematic segregation of galaxies with 
diverse specific-SFRs is apparent. Galaxies with the highest 
specifle-SFRs (>10“® ®yr“'; blue cross) have the highest ve¬ 
locity dispersion as a population, and a higher mean cluster¬ 
centric radius, than those star-forming galaxies with the low¬ 
est specific-SFRs (< 10“'° ^ yr“'; red cross). There is a general 
progression towards lower mean radii and LOS velocity dis¬ 
persions with decreasing specific-SFR. This progression con¬ 
tinues to the inactive galaxy cluster galaxy population (not de¬ 
tected at 24/2m; black cross), which has a significantly lower 
mean radius and velocity dispersion than any of the four sub¬ 
populations of star-forming galaxies. 

The mean stellar masses of the star-forming cluster galax¬ 
ies increases by 0.5 dex from the highest specific-SFR bin to 
those in the lowest one, reflecting the systematic decline in 
spe cific-SFR with stell ar mass for star-forming cluster galax¬ 
ies (iHaines et al.ll2013T) . We may thus be concerned that these 
kinematical differences are in fact due to mass segregation 
rather than a sequence in declining specific-SFR. 

There is no evidence for mass segregation within 
our cluster galaxy popul ation however, consistent with 
Ivon der Linden et al.l (120101) . Splitting the cluster galaxy pop- 



mean cluster-centric radius < r > / rjoo 

Figure 8. Variation in the spatial distribution and kinematics of cluster galax¬ 
ies as a function of their specific-SFR and stellar mass. Colored crosses 
indicate the mean cluster-centric radii {{rpwj/fsoo)) and LOS velocity dis¬ 
persions {cr/o.y(rp^o;< 2 . 0 r 5 oo)/<^i^) of star-forming cluster galaxies split into 
four bins of specific-SFR. En'or ellipses show the Icr uncertainties in both 
values. The black cross indicates the corresponding values for inactive clus¬ 
ter galaxies not detected with Spitzer. Only cluster galaxies within 2r5oo are 
included within each sub-sample. The mean stellai* mass, {log(vM/MQ)), 
of each sub-sample is indicated in the lower-right comer. Thin dashed el¬ 
lipses indicate the results after applying additional stellai* mass cuts to each 
specific-SFR bin to equalize their mean stellar masses to that of the inac¬ 
tive population. Grayscale symbols and error-bars indicate the mean cluster¬ 
centric radii and LOS velocity dispersions of cluster galaxies (including both 
star-forming and inactive sub-populations) split into four bins of stellar mass: 
log(Ad/M0)<lO.4 (light gray circle); 10.4<logA4 <10.6 (gray triangle); 
10.6< log 10.85 (mid-gray diamond); log Ad >10.85 (dark-gray square). 

ulation into bins of stellar mass {grayscale symbols), much 
less variation in the mean radii and LOS velocity dispersions 
is seen between stellar mass bins in comparison to those split 
by specific-SFR, and no overall trend with stellar mass is vis¬ 
ible. Moreover, the kinematic segregation by specific-SFR 
persists even if additional stellar mass cuts are applied to each 
specific-SFR bin to equalize their mean stellar masses {thin 
dashed ellipses). 

4. MAPPING THE CONTINUAL ACCRETION OF GALAXIES ONTO 
MASSIVE CLUSTERS IN THE MILLENNIUM SIMULATION 

To correctly interpret the previous observed trends in clus¬ 
ter galaxy properties, clusters must be placed in the cosmo¬ 
logical context of continually accreting galaxies and groups 
from their surroundings. With this aim, we have examined 
the spatial distributions and orbits of galaxies in the vicinity 
of th e 75 most massive cl usters from the Millennium simula¬ 
tion (ISpringel et al.ll2005h . a cosmological dark matter simu¬ 
lation covering a (500/i“'Mpc)^ volume. These clusters have 
present day virial masses in the range 4.0-23.6x 10'"'/i“' Mq, 
velocity dispersions of 630-1540 km s“', and a median forma¬ 
tion epoch, Zf=0.59. We have extracted 20x20x 140/2“^Mpc^ 
volumes centered on each cluster. These volumes are ex¬ 
tended in the z-direction so that, for a distant observer view¬ 
ing along this axis, all galaxies with LOS velocities within 
5 000 km s“' of the cluster redshift are included, enabling pro¬ 
jection effects to be fully account for and quantified. 

There exists a full database of properties for each galaxy 
in the simulation, including positions, peculiar velocities, ab- 
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Redshift at which galaxy is accreted by cluster 

Figure 9. The orbits of galaxies about the 10th most massive cluster in the 
Millennium simulation, as viewed by a distant observer along the z-axis. The 
orbit of each galaxy with A1>2xiO^^MQ at z=0 is shown by a colored 
curve, tracing its movement from z=0.76 (snapshot 44) to z=0.0 (snapshot 
63). The final location of the galaxy is marked by a dot, while the cross indi¬ 
cates its location at z=0.41 (snapshot 50). Each curve is color coded accord¬ 
ing to the epoch at which the galaxy is accreted into the cluster, as indicated 
at the bottom of the plot. This accretion epoch is defined as the snapshot at 
which the galaxy passes within ^ 200 ( 2 ) for the first time. Galaxies yet to pass 
within r 2 oo have mid-blue colors, while those accreted earliest into the clus¬ 
ter (Zflcc>0.51) have red colors. The dashed black circle indicates the present 
day ^200 radius of the cluster. Only galaxies which would be spectroscopi¬ 
cally identified as cluster members by the observer are shown. 

solute magnitudes, stellar masses etc., b ased upon th e GAL- 
FORM semi-analytic models (SAMs) of iBower et al.l (120061) 
at 63 snapshots throughout the life-time of the Universe to 
z=0. Similarly, another database provides the positions, veloc¬ 
ities, masses (M 200 ) and radii (r 2 oo) of each DM halo at each 
snapshot. For each galaxy and halo in a given snapshot, the 
database provides links to identify its most massive progenitor 
in the preceding snapshot, and so on, all the way back to its 
formation, allowing its mass growth and full merger history to 
be mapped in detail (see e.g. De Lucia & Blaizot 2007). This 
process also allows the orbit of each galaxy with respect to 
the cluster to be followed from formation to the present day, 
enabling us to determine its epoch of accretion (Zacc) into the 
cluster, defined here as the redshift at which the galaxy passes 
within r 2 oo(z) for the first time. 

4.1. The infall of galaxies onto clusters 
Figure |9] shows the orbits of galaxies about the tenth most 
massive cluster (M2oo=9.9x 10*"^/i“' Mq at z=0) in the Millen¬ 
nium simulation, color coded according to accretion epoch. 
Almost all the galaxies over the 20x20/i“^Mpc^ field-of-view 
are falling steadily into the cluster or have already been ac¬ 
creted. The complex large-scale structure around the cluster is 
apparent, including clear preferential directions for the galax¬ 
ies to flow into the cluster, while other regions appear largely 



Look-back time (Gyr) 


Figure 10. The infall of galaxies onto a massive cluster. Each galaxy with 
Ad>2xIO'^M q at z=0 is shown by a colored curve tracing its comoving 
cluster-centric distance as a function of look-bacic time. The curves are color 
coded according to the epoch at which the galaxy is accreted into the clus¬ 
ter, as in Eigure|2] The first pericenter and apocenter of each galaxy’s orbit 
about the cluster are marked with dots. The black dashed curve indicates the 
evolution of the cluster radius r 2 oo with time. Only galaxies which would be 
spectroscopically identified as cluster members by the observer are shown. 

devoid of galaxies. Galaxies tend to be drawn first into the 
filaments from the surrounding field, then flow along the fil¬ 
aments into the cluster. While many galaxies are infalling as 
individual objects from the field, others are arriving into the 
cluster within galaxy groups (e.g. the tangle of green curves 
coming in from the right-hand side). 

The full extent of the cluster’s gravitational sphere of influ¬ 
ence is revealed in Fig.[T0] with all galaxies within 18/r*Mpc 
of the cluster falling steadily inwards. For all 75 clusters, the 
boundary between the infall regions and beyond, where galax¬ 
ies remain attached to the Hubble flow, is found at a comoving 
distance ^10-20/i“'Mpc from the cluster. 

The infall of galaxies into the cluster is highly coherent, at 
least for z^l: the radial velocities of infalling galaxies at the 
same cluster-centric radius at a given epoch are all roughly the 
same, and the future trajectories and accretion epochs of an in¬ 
falling galaxy can be accurately estimated simply on the basis 
of their current cluster-centric distance. For example, almost 
all galaxies which were 5/!“'Mpc from the cluster 4 billion 
years ago, are due to be accreted into the cluster at z^O.l, as 
indicated by the parallel diagonal green colored curves. 

After accretion, galaxies remain bound to the cluster, but 
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Figure 11. Radial phase-space diagram (Wradial vs r/rioo) for galaxies with 
A4>5 X10^ Mq orbiting around the same massive cluster as in Figsl9land ll()l 
Galaxies are color coded according to their accretion epoch as indicated. 

many have orbits which take them outside r 2 oo {black dashed 
curve), including some which bounce out as far as ~2.5r2oo. 
These galaxies, rebounding out of the cluster after their first 
pericenter passage are known as the “back-splash” popula¬ 
tion. The time-scale required for galaxies to reach pericen¬ 
ter after being accreted is of the order 0.5-0.8Gyr, while 
the orbital periods of the “back-splash” galaxies are much 
longer, only reaching their first apocenter 2-3 Gyr after pass¬ 
ing through the cluster for the first time. 

4.2. Phase-space diagrams 

Figure[TT]shows the distribution of galaxies in radial phase- 
space: radial velocity {vradud) versus cluster-centric radius (r) 
in units of raoo, for the same cluster, where Vradiai is the ra¬ 
dial component of the galaxy’s velocity relative to the cluster, 
including a component from the Hubble expansion {Hq r). 

The distribution of galaxies in phase-space splits into two 
reasonably well defined structures: a triangular-shaped virial- 
ized region containing galaxies which have passed through 
the cluster at least once; and a narrow stream of infalling 
galaxies {blue p oints) with negative radial velocities ext end- 
ing out to 4 r 2 oo (iMamon et al. [2004 iDiinner et al.ll200'^ . As 
the infalling galaxies plunge into the deep gravitational po¬ 
tential well of the cluster, they are continually accelerated, 
reaching infall velocities of up to 3 000 km s“'. After passing 
through pericenter, these galaxies reappear along the top edge 
of the triangular region as back-splash galaxies, coherently 
progressing outwards and slowing down with increasing Zacc 
towards the right-hand apex at >2 r 2 oo that marks the apocen¬ 
ter of the orbits of those galaxies accreted ~3 Gyr ago. The 
radial phase-space diagram retains much of the information 
regarding the epoch of accretion of a galaxy, allowing this 
epoch to be accurately estimated for galaxies based on their 
location in the diagram, at least for those accreted within the 
last 3 Gyr. Only those which were accreted much earlier have 
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Figure 12. The caustic diagram (Aiz/o, vs rpr„jlr 2 oo) for M>2x 10*®Mq 
galaxies orbiting around the same massive cluster, as viewed by distant ob¬ 
server along the z-axis. Galaxies are color coded according to their accretion 
epoch as in Figs. njtn] Galaxies within 20 h ’ Mpc of the cluster center at 
z=0.0 are indicated by solid symbols, while those at comoving cluster-centric 
radii r>20/i“* Mpc are shown by open symbols. The velocity dispersion of 
galaxies within a projected cluster-centric radius rproj<r 2 Qa is indicated. 

had time for their orbits to become mixed. 

Figure [T2| shows the observable counterpart to Fig.[TT] the 
caustic diagram, which plots the LOS velocity of galaxies rel¬ 
ative to the cluster redshift (Atz/oj) against projected cluster¬ 
centric radius {rproj/r 2 m) for the same cluster, as viewed by 
a distant observer along the z-axis. The relative LOS velocity 
of galaxies combines the LOS component of their peculiar ve¬ 
locities with the contribution to their redshifts from the Hub¬ 
ble expansion: Avios={i^z,gai-v^,ci) +Hodz, where Uz,gai and 
idz,cl are the LOS peculiar velocities of the galaxy and clus¬ 
ter DM halo respectively, and dz is the distance between the 
galaxy and the cluster halo along the line of sight. 

The distribution of cluster galaxies shows the typical 
“trumpet”-shaped caustic profile. Throughout this paper we 
refer to “spectroscopic” cluster members as those galaxies 
which lie within the caustic profiles determined for each of 
the 75 clusters. This is not the same as galaxies which are 
identified as satellites of the cluster DM halo, which are those 
physically located within the cluster halo (r<r 2 oo), or even 
the population of accreted galaxies that have passed within 
f 20 o{z) at some point in their history (but which may now be 
outside r 2 oo). For this cluster, the separation between those 
galaxies that would be identified as “spectroscopic” cluster 
members, and those which are clearly background objects, 
roughly corresponds to a separation between those physically 
within 20/1“* Mpc of the cluster center {solid points), and 
those beyond this radius {open points). Indeed, all the galaxies 
shown earlier in Figure [T0| are “spectroscopic” cluster mem¬ 
bers within a projected cluster-centric radius 2.0r2oo- Exam¬ 
ining each of the 75 clusters individually, the cluster-centric 
distances that best separate “spectroscopic” members and 
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Figure 13. Stacked phase-space diagram (ctu vs. rproj /^ 200 ) for the 75 most massive clusters in the Millennium simulation at z=0.0 (panel a). Each point 
indicates an Ad >2 x 10^^ M 0 galaxy from the lBower et all 420Q6D semi-analytic model catalog, colored according to when it was accreted onto the cluster. The 
four panels on the right split these cluster galaxies into dynamical sub-populations: (b) just those galaxies yet to pass within r 2 oo and be accreted into the cluster; 
(c) infalling galaxies {i'radiai<^) which have passed within r 2 oo^ but have yet to reach pericenter; (e) back-splash galaxies which have passed through pericenter 
and are on their way back out of the cluster {Uradiai >0) and are yet to reach apocenter; (f) the virialized cluster population which have passed through apocenter. 
The lower left panel (d) displays the velocity dispersion profiles {o'ios{f')/cFij) of each sub-population as well as that of the overall cluster population. 


clear fore/background objects are 19±3/i“* Mpc (13±2r2oo)- 

Unlike the radial phase-space diagram, it is not possible to 
select individual galaxies from a specific region of the caustic 
diagram and then identify it as an infalling, recently accreted 
or virialized galaxy. However, several trends within the dis¬ 
tribution of galaxies in the caustic diagram can be seen. First, 
those galaxies accreted earliest (Zace^O.4; red points) are spa¬ 
tially localized in the cluster core with typical LOS veloci¬ 
ties <1 000kms“*. Second, beyond rpr„j>l.S r 2 oo the bulk of 
galaxies have yet to be accreted into the cluster while the re¬ 
mainder all appear to be back-splash galaxies accreted ~3 Gyr 
ago. Finally, many of the galaxies with the largest LOS ve¬ 
locities (> 1 OOOkms”') have only recently arrived into the 
cluster (green points). Although these general trends hold for 
the vast majority of the clusters in our sample, there is signifi¬ 
cant cluster-to-cluster variation in the spatial distributions and 
relative contributions of galaxies accreted at different epochs, 
as expected given their dynamical immaturity. 

4.3. Stacking the clusters 

To account for the cluster-to-cluster scatter in a statistical 
way, the caustic diagrams for all 75 clusters are stacked to 
produce Fig. [T3h . The cluster-centric radius of each cluster 
member is scaled by the rioo of that cluster, and the LOS ve¬ 
locities are scaled in units of crv(rproj<riao), the LOS velocity 
dispersion of all “spectroscopic” members within r 2 oo- 


To demonstrate how the distribution of galaxy populations 
in the caustic plot provides information about their dynamical 
evolution and accretion history, the right-hand panels show 
the caustic diagrams after splitting the cluster population into 
four dynamical sub-populations, while panel (d) plots each 
of their velocity dispersion profiles (Tios(r)l<Ji, alongside that 
for the overall cluster population. First, galaxies yet to pass 
within r 2 oo and be accreted (panel b) are found at all cluster¬ 
centric radii, becoming increasingly dominant with radius. 
The fall in numbers towards the cluster core (rproy^O.b r 2 oo) is 
due to the area of sky covered in any given narrow radial slice 
scaling as r, rather than a physical decline in the surface den¬ 
sity of such galaxies at low radii (see Fig. [nil. The velocity 
dispersion of these galaxies remains relatively constant with 
radius, as they have yet to approach the cluster core 

Panel (c) shows those infalling galaxies which have recently 
passed within r 2 oo for the first time, but have yet to reach 
pericenter. These objects are all found within r 2 oo by de¬ 
fault. They have the highest LOS velocities of all our sub¬ 
populations (cr/os(r)> 1.4 cTjy), being fully accelerated as they 
fall deep into the gravitational potential well of the cluster 
core, and at late epochs when the clusters are much more mas¬ 
sive than at any previous point in their history. 

Panel (e) shows the back-splash population of cluster galax¬ 
ies which have passed through pericenter, and are now head¬ 
ing back out away from the cluster center ii'radiai>0), but have 
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Figure 14. Radial population gradients in clusters. Each cui've shows the fraction of “spectroscopic” cluster galaxies which have been accreted by a given 
redshift as a function of projected cluster-centric radius (normalized by ^ 200 ) averaged over our stacked sample of 75 massive clusters observed at z=0.0, color 
coded according to accretion redshift as indicated on the far right of the plot. The three panels indicate the trends produced when considering three definitions 
for describing when a galaxy is accreted into the cluster, taken to be the snapshot when the galaxy passes within rsoo panel), ?'200 {middle panel) or 2?'200 
{right panel). The vertical dashed line in each panel indicates the corresponding radius at which the accretion epoch for the curves in that panel are defined. 


yet to reach apocenter. These galaxies were typically accreted 
l^Gyr ago. This population shows a triangular distribution 
in the caustic diagram, with high LOS-velocities (tr/o^^LS- 
1.5 (Tj/) in the cluster core, having recently completed their 
first infall, which steadily fall as the galaxies rebound out 
of the cluster, slowing down as they attempt to climb back 
out of its potential well. The back-splash population ex¬ 
tends as far out as 2-3 r 2 oo, where they can be differentiated 
from the infalling population by their characteristically low 
LOS-velocity dispersions (cr/o 5 (r)^ 0 . 4 - 0 . 55 cr,y), some ~35% 
lower than the overal l cluster population at the sa me cluster¬ 
centric distances (see lGill. Knebe & Gibsonl[2005l) . 

Finally, the virialized population (panel f) are those galaxies 
which were accreted at early epochs (Zacc^O.4) and have all 
passed through apocenter of their first orbit. This population 
also presents a triangular distribution in the caustic diagram, 
but is much more concentrated towards the cluster core (par¬ 
ticularly at < 0 . 2 r 2 oo) than the back-splash galaxies, and has 
lower LOS-velocity dispersions at all radii. We identify this 
population with those galaxies that either formed locally or 
were accreted when the cluster’s core was being assembled. 
Their low velocities reflect the fact that the system they fell 
into was much less massive than the present-day cluster. 

4.4. Radial population gradients 

Figure [14] shows the radial population gradients obtained 
when stacking the 75 massive clusters as observed at z=0.0. 
The thick blue curve in the middle panel shows how the frac¬ 
tion of “spectroscopic” cluster members which have been ac¬ 
creted varies as a function of projected cluster-centric radius, 
rproj / 1 ‘ 200 . This fraction drops slowly from 100% at the cluster 
core to ^90% at 0 . 4 r 2 oo, before falling at an ever increasing 
rate down to ^50% by r 2 oo, and approaching zero by ^3 r 2 oo. 
For rproj>r 2 oo, this contribution represents the “back-splash” 
population galaxies which have previously been accreted into 
the cluster, but now have bounced back out beyond r 2 oo. 

The remaining curves show the effects on this fraction by 
progressively excluding galaxies which were accreted after a 
given redshift. The first three curves (zacc<0.09) only diverge 
from the top curve at rproj <r 2 oo, as the most recently accreted 
galaxies haven’t had sufficient time to pass though the clus¬ 
ter and go back out beyond r 2 oo. As those galaxies accreted 
1-3 Gyr ago are progressively excluded, the curves become 
increasingly steep within r 2 oo, and flatter beyond r 2 oo. The 


bulk of the “back-splash” population found beyond r 2 oo were 
accreted between z=0.14 and z=0.51. 

The remaining panels show the effect of changing the nom¬ 
inal radius for identifying when a galaxy has been accreted 
into a cluster. Reducing the accretion radius from r 2 oo to rsoo 
(left panel) simply squeezes the curves inwards by a compa¬ 
rable amount. Pushing the accretion radius outwards to 2 r 2 oo 
{right panel) increases the prominence of the bump produced 
by recently accreted galaxies, but greatly reduces the contri¬ 
bution from back-splash galaxies, as so few galaxies which 
pass within 2 r 2 oo rebound from the cluster beyond this radius. 

The key aspect of all these curves is that irrespective of ac¬ 
cretion epoch, Zacc, or the precise radius used to define ac¬ 
cretion (within reasonable limits), the resulting radial popu¬ 
lation gradient drops from ^ 100 % in the cluster core to zero 
by ~ 3 r 2 oo. This steep gradient is primarily due to the com¬ 
plementary increase with radius in the fraction of “interloper” 
galaxies that have yet to be accreted into the cluster among 
the “spectroscopic” cluster population, from zero in the clus¬ 
ter core to 100 % by ~ 3 r 2 oo. A second contribution to the 
steepness comes from the correlations between the epoch of 
accretion of satellite galaxies a nd their present cluster-centric 
distance (iDe Lucia et akl 12013) . with those satellite galaxies 
close to the cluster core having been accreted significantly 
earlier on average than those located close to the viria l ra¬ 
dius (both physically or in projection). iGao et all (l2004ft find 
the same radial trend for sub-halos, with the median accretion 
redshift of sub-halos in massive cluster halos decreasing from 
Zacc^^-0 in the cluster core to Zacc^OA at r 2 oo (their Fig. 15). 

5. CONSTRAINING STAR FORMATION QUENCHING MODELS BY 
COMPARISON TO THE OBSERVED SF-RADIUS TRENDS 

Figures |2]BI showed that the fraction of massive cluster 
galaxies with SFR/«> 2 .OM 0 yr“' or blue UV-optical colors 
increases steadily with cluster-centric radius, but at the largest 
radii probed, the fsp remained significantly (20-30%) below 
that seen in coeval field galaxies in both cases. 

In Figure [15] we attempt to reproduce these two SF-radius 
trends by a simple toy model in which the star formation of 
infalling field galaxies is instantaneously quenched at the mo¬ 
ment they pass within r 2 oo of the cluster for the first time, or 
after a certain time delay (At). The fraction of star-forming 
galaxies among this infalling field population is set to match 
our observed coeval field galaxy sample. The stacked radial 
population gradients for Af> 2 xlO'®M 0 galaxies from the 
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Figure 15. Comparison of the observed radial population gradients with 
predictions from cosmological simulations. The red points show the frac¬ 
tions of Ad>2xlO*^M0 cluster galaxies with obscured SFR>2.0 Mq yr“^ 
from Fig. [21 while the magenta squares show the corresponding SF-radius 
relation for UV-selected star-forming galaxies (Fig. |4j. The blue horizontal 
line indicates the corresponding fraction of field galaxies (Ad>2x IO^^Mq; 
0.15<^<0.30) with SFR/^ > 2.0MQyr“^ while the blue shaded region in¬ 
dicates the Icr confidence levels. The thick blue diagonal curve shows the 
predicted SF-radius relation obtained from our “observations” of the 75 mas¬ 
sive clusters in the Millennium simulation, assuming that infalling galaxies 
have the same fsF as our observed field galaxy sample, and star-formation 
is immediately quenched upon being accreted into the cluster, i.e. passing 
within ^200 for tho first time. The colored curves show the effects of delaying 
this quenching by a time At after the galaxy is accreted into the cluster, as 
indicated by the color scale below the plot. The dashed blue curve shows 
the effect on the predicted SF-radius relation of changing the cluster-centric 
radius at which quenching is initiated from r 2 oo out to 2 r 2 oo- 

same 75 massive clusters in the Millennium simulation are 
reproduced, as they would appear if observed at z=0.21, to 
best match the redshifts of the LoCuSS sample. The model 
galaxy positions and velocities relative to the cluster halo are 
now measured as they stood at z=0.21, while the clusters are 
stacked using their r^oo and ajy values measured at z=0.21. 
At z=0.21 these clusters have M 200 masses in the range 2.6- 
21.7x 10*4/r' Mq, with a median Mjoo of 5.Ox 10' V' Mg. 

The predicted SF-radius relation in the case that star for¬ 
mation in all infalling galaxies is instantaneously quenched 
upon accretion {thick blue diagonal curve) is qualitatively 
similar in form to the observed trends, and consistent with 
the data points at ~l- 2 r 5 oo, suggesting that this is to 
first order a reasonable as s umptio n, as found previously b y 
iBalogh. Navarro & MorrisI (l2000h and iHaines et~aD (l2009h . 
The model radial gradient is too steep however, resulting in 
predicted values of fsF that are much higher than our data 
points in the range 1.8-3.0 rsoo, and too low in the cluster core 
(VofZO-8/‘50o)- 

The remaining colored curves show the effects of delay¬ 
ing the moment at which quenching occurs, by terminating 
star-formation only in those galaxies accreted into the clus¬ 
ter more than Af Gyr prior to observation, corr esponding to 
the “delayed-then-rapid” quenching scenario of IWetzel et alJ 
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Figure 16. Radial population gradients. The red and magenta poi nts show 
the observed IR- and UV-based SF-radius r elat ions as in Figure flSl The col¬ 
ored curves are the same as those in Figure 1151 except that they now assume 
that the infalling galaxies have a fraction fsF which has been reduced by 19% 
with respect to that seen in the general field (blue horizontal line). The black 
dot-dashed curve indicates the predicted SF-radius relation in the case that 
star formation is instantaneously quenched when galaxies reach pericenter. 

dMl- The ‘ ‘excess” obscured star-formation {redpoints) ob¬ 
served in the cluster core can then be reasonably reproduced 
by a model with a short quenching delay of the order 0.3- 
1.0 Gyr {light-blue/green curves). Much longer quenching 
time-delays (Af>3 Gyr) are clearly excluded, as they leave 
too many star-forming galaxies at - 2 r 5 oo- 

One possible way to reconcile the model predictions with 
our data at ^3-5 7500 would be to initiate the quenching pro¬ 
cess at larger radii. The dashed curve shows the radial popu¬ 
lation gradient produced when star-formation in all infalling 
galaxies is immediately quenched when they pass within 
27200 for the first time (rather than 7200 as before). While the 
fraction of star-forming galaxies has now been reduced suf¬ 
ficiently at ^ 27500 , the gradient of the model trend remains 
much steeper than that observed, and vastly under-predicts the 
fraction of star-forming galaxies at rproj^l.5 7500 . 

The model curves provide little or no leeway to reproduce 
the low fsF observed over 3-5 7500 in the NUV-based SF- 
radius relation {magenta squares), as they all approach the 
field value at 3.5—4.0 7500 ("^2.5-3.072oo). the radius at which 
the fraction of back-splash galaxies falls to zero (Fig.fT4li. ir¬ 
respective of the radius at which quenching is initiated. 

The only way to improve the model fits to these observa¬ 
tions would be to allow the fraction of star-forming galax¬ 
ies among the infalling population to be lower than that ob¬ 
served among the general coeval field population. We con¬ 
sider the simplest approach in Figure [16] that of simply re¬ 
ducing the fsF of cluster galaxies by a single fixed amount 
(19%) from the value observed in the field, at all radii, to 
model the impact of whatever physical process is reducing 
star formation among galaxies in the infall regions of clus¬ 
ters. One feasible mechanism to achieve this is via the “pre- 
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processing” of galaxies within galaxy groups which are subse- 
q uently accreted into clu sters. The cosmological simulations 
of lGabor & Dav^ (12014^ suggest that ^40% of satellite galax¬ 
ies within clusters are pre-processed, the fraction decreasing 
weakly with cluster-centric radius. The overall IR-based SF- 
radius relation can now be reproduced at all radii, via a model 
in which star-formation is quenched in galaxies ^0.7-1.5 Gyr 
after being accreted into the cluster. This occurs on average 
after the galaxy has passed through the cluster core, as the ob¬ 
served trend lies above that predicted for the case in which star 
formation is quenched in galaxies at the moment they reach 
pericenter (black dot-dashed curve). 

The NUV-based SF-radius relation is systematically above 
that of the IR-based inside ~2.5 rsoo, and is best-fit by a model 
is which star-formation is quenched in galaxies ^2.1-3.6 Gyr 
after accretion, or slightly later than that suggested by the IR- 
based relation. 

For all t hese model curves, to use the terminology of 
iPeng et all (l2()Tnl [2?iTll . we assume that the environmental 
quenching process is 100 % efficient (eja,=l), i.e. all star¬ 
forming galaxies are quenched Af Gyr after accretion into the 
cluster. The result of decreasing Csat for clusters would leave 
the model curves unchanged at large radii (ry„oj> 3 r 2 oo), but 
squeeze the curves upwards in the core, effectively reducing 
the radial population gradient by a factor (1 - Csar)- The fact 
that the observed gradients are steep, with fsF increasing five¬ 
fold over 0-3 r 2 oo, necessitates a high environmental quench¬ 
ing efficiency, e 5 a,> 0 . 8 , otherwise there would be evidence of 
a residual population of primordial cluster galaxies with on¬ 
going star-formation in cluster cores (rpmj^O-^ i' 20 o)- 

5.1. Uncertainties in the model trends 

One concern about attempting to fit the SF-radius relation 
with model curves, is that the latter require the semi-analytic 
models to reliably follow the orbits of cluster galaxies long 
after their accretion, and accurately model the long-term evo¬ 
lution of their stellar masses once they become satellites. The 
Millennium simulation itself considers only the dark matter 
component, the baryons bolted on afterwards via SAMs. 

When galaxies are accreted into massive clusters, their par¬ 
ent DM halos become sub-halos of the cluster halo, losing 
mass continuously through tidal stripping, until in many cases 
they fall below the resolution limit of the simulatio n, dissolv¬ 
ing into the parent halo or are completely disrupted (iGao et alJ 
l2004t IWeintnann et al.ll2010h . The galaxies hosted by these 
sub-halos are much more compact and tightly bound than the 
dark matter, and as long as some of the surrounding sub¬ 
halo survives, are not expected to suffer significant stellar 
mass loss, although their diffuse hot gas halo is likely to be 
lost at a rate com mensurate with that of the parent sub-halo 
(iGuo et al.ll201 ll) . Prior to their accretion into the clusters, 
model galaxies at our lower stellar mass limit of 2 x IO'^Mq 
typically have parent DM halos of masses ^ 6 x lO'*/r'M 0 in 
the Millennium simulation, comprising ~730 DM particles, 
and hence must suffer > 95% stripping before their parent 
sub-halo falls below the mass resolution limit of 20 particles 
( Snringel et al.ll2005[). a process which typically takes ^5 Gyr 
( Weinmann et al.ll2010t) . Even after their parent sub-halo has 
been entirely disrupted, they are expected to survive as “or¬ 
phan” galaxies, although now they will be likely subject to 
significant stellar mass loss via tidal stripping and may be 
completely disrupted. 

The inclusion and treatment of “orphan” galaxies by SAMs 
is required to explain the cuspy radial satellite number den¬ 


sity profiles observed within cluster halos (Fig. O, which 
are much steeper than the relativel y flat radial distribu- 
tions of DM sub-halos within 0.3 rsoo (IBudzvnski et al.ll2012l : 
IVogelsberger et al1l2014f) . These SAMs resort to ad hoc pre¬ 
scriptions for mass-stripping and adjustment of the orbits of 
these “orphan” galaxies within the cluster halos, resulting in 
strong variations in t he radial satellite numb er density profiles 
within cluster cores (IBudzvnski et al.ll2012h . and the fractions 
of “orphans” i n the cluster satellite popu lation, from 25% in 
the lGuo e t al.l(l201 ll) model to 50% in the lBower et al.l (120061) 
model dGifford etffll2013l) . This likely explains the inabil¬ 
ity of th e predicted radial p rofile of model cluster galaxies 
from the iBower et al.l (I2006h SAM to match the observed ra¬ 
dial profile at < 0.1 rsoo, overestimating the number den¬ 
sity of cluster galaxies by a factor 2-3 (Fig.|5]l. 

We should therefore be cautious about using model curves, 
such as those in Figs. [15] and [16] that depend upon the radial 
distribution and numbers of those galaxies accreted earliest 
into the clusters, and which have suffered repeated interac¬ 
tions over multiple orbits within the ever growing cluster halo. 
All the model curves in these two figures do this, even those 
referring to galaxies yet to be accreted into cluster, as they de¬ 
pend upon the relative contributions of those accreted after a 
given epoch, with all those accreted before the same epoch. 

5.2. Radial galaxy surface density profiles 

Figure [it] compares the radial galaxy surface density pro¬ 
files, E(r), of star-forming galaxies with model surface den¬ 
sity profiles considering just those galaxies accreted within 
the last Af Gyr or which have not yet been accreted. This plot 
largely resolves the above issues, by focusing solely on star¬ 
forming galaxies, which are most likely to still have surviving 
parent sub-halos, and comparison model curves that contain 
just the most recent arrivals into clusters, and hence minimiz¬ 
ing the uncertain contribution from “orphan” galaxies. 

The surface density of star-forming (SFR/s> 2 .OM 0 yr“*) 
cluster galaxies from our ensemble of 30 clusters (blue points) 
declines steadily with radius out to ^3 rsoo, with no evidence 
of flattening off inside rsoo. This immediately rules out mod¬ 
els in which star formation is instantaneously quenched when 
galaxies are accreted into clusters (thick blue curve), as these 
produce radial profiles which are essentially flat within 2 rsoo- 
As recently accreted galaxies are progressively included, the 
radial density profile steadily builds up and steepens within 
2 rsoo- The best fitting curve to the Spitzer data is obtained 
by considering those model cluster galaxies accreted within 
the last Af=2.1j;Q 7 Gyr, with a value of 21.15 for 23 data 
points and two degrees of freedom (f and Af). The uncertain¬ 
ties in Af are derived as the values for which the value has 
increased by 2.30 from the minimum value. 

Outside of 0.3 rsoo, the shape of the Y,(r) profile 
for unobscured star-forming galaxies (magenta squares', 
Al> 2 xlO*°M 0 , (NUV - r)o,Q < 4.5), coincides well with 
that of the obscured star-forming population. In contrast, the 
E(r) profile for the UV-selected star-forming galaxies steep¬ 
ens more rapidly in the cluster core to form a cusp, parallel¬ 
ing that seen for the overall cluster population (red points). 
The best fitting curve to the NUV data has Af=3.2±0.4Gyr 
(X^=29.57 for 32 data points), implying that the NUV emis¬ 
sion takes 1.1 Gyr longer after accretion to be shut down than 
the 24p,m emission. 

Figure [T 8 ]replots the E(r) profile for obscured star-forming 
cluster galaxies (blue points), but now compares it with slow 
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Figure 17. Blue symbols show the S(r) profile of A4>2x IO'^Mq galax¬ 
ies with obscured star formation at rates SFR//;>2.OM0yr“* averaged over 
our 30 clusters. Magenta squares show the S(r) profile for unobscured star- 
forming galaxies with (NUV—r)o,o<4.5, normalized to match the profile of 
obscured star-forming galaxies at rprDj>l -5 rjoo. The thick blue curve shows 
the predicted surface density profile of infalling galaxies yet to pass within 
floo, obtained from our “observations” of the 75 massive clusters in the Mil¬ 
lennium simulation at z=0.2L noimalized to fit the observed radial profiles of 
star-forming galaxies at large radii. The remaining colored curves show the 
predicted radial profiles produced by including also those galaxies accreted 
into the clusters within the last At Gyr as indicated by the color scale. The 
black curve shows the predicted surface density profile for all “spectroscopic” 
member galaxies of the same 75 clusters, scaled to best match the observed 
ensemble S(r) profile of all Ad>2x IO'^Mq cluster galaxies (redpoints). 

quenching models. In this scenario model star-forming galax¬ 
ies are accreted into the cluster and subsequently gradually 
quenched, their SFRs declining exponentially on a quench¬ 
ing time-scale fg until their SFRs fall below our nominal 
limit of 2.0M0 yr“*. These model star-forming galaxies are 
given initial SFRs taken at random from our observed sam¬ 
ple of coeval field star-forming galaxies (SFR>2M0yr“'; 
0.15<z< 0.30), the SFR distr ibution of which is shown in 
Fig. 2 of iHaines et al.l (1201 3 ll . This should be reasonable 
given that the infrared luminosi ty functions of cluster and field 
galaxie s are in distinguishable (iFinn et al.ll2010t IHaines et al.l 
I2OI lallR I2OI3I) . and the specific-SFRs of star-forming galax¬ 
ies in infall re gions are indistingu ishable from those in coeval 
field samples iHaines et al.ll2013]) . The colored curves show 
the predicted surface density profiles for quenching time- 
scales fg in the range 0.25-10.OGyr. The best fit model to 
observations has a quenching time-scale fQ=2.19±0.41 Gyr, 
with a value of 20.80 for 23 data points and two degrees of 
freedom (</>, tg). 

5.3. Dynamical analysis 

The velocity dispersion profile (VDP) of cluster galaxies, 
o'losir), provides complementary constraints for the accretion 
epochs of galaxy sub-populations (see Fig. fT3ll). Figured^ 


Figure 18. Blue symbols show the S(r) profile of Af >2x lO’® M0 galaxies 
with obscured star formation at rates SFR;x'>2.OM0yr“*. Red symbols in¬ 
dicate the con'esponding S(f') profile for all Ad >2 X10*® M© galaxies. The 
thick blue curve shows the predicted surface density profile of infalling galax¬ 
ies yet to pass within r 200 . normalized to fit the observed surface density 
profile of star-forming galaxies at large radii. The remaining colored curves 
show the predicted surface density profiles for star-forming galaxies accreted 
into the clusters and subsequently have their star-formation rates decline ex¬ 
ponentially with a range of quenching time-scales (rg=0.25-10 Gyr). 

compares the observed VDPs of 24/rm-detected star-forming 
galaxies {blue points) with the predicted VDPs of model clus¬ 
ter galaxies, selected according to their accretion epoch. 

The VDP for star-forming cluster galaxies shows a high, 
narrow peak of 1.44cry at Pproj ^ 0.3 7500, before dropping 
to the innermost radial bin, and a steady decline outwards to 
^0.8 (Ti/ at large radii. This profile shape is best reproduced 
by model cluster populations combining infalling galaxies yet 
to be accreted and the most recent arrivals into the cluster. 
The progressive inclusion of these recently accreted galax¬ 
ies causes the velocity dispersion within 7590 to rise rapidly, 
producing a characteristic sharp peak at 7^^7^0.27500 which 
reaches a maximum height of ^1.56cri/ when At = 0.7 Gyr, 
along with a corresponding sharp drop off to the cluster core. 
The best overall match to observations is produced for mod¬ 
els with Af^0.5-2.2 Gyr {light blue/green curves), compara¬ 
ble with the time-scales required for infalling galaxies to ap¬ 
proach the pericenter of their orbits through the cluster and 
achieve the high velocities required to produce the observed 
peak in <Tios{r). The observed profile is inconsistent with 
models with much longer delays between accretion into the 
cluster and quenching (Af>3.7Gyr), and models in which 
star-formation is quenched instantaneously upon accretion 
(Af=0 Gyr; thick blue curve), due to their predicted low, rela¬ 
tively flat LOS VDPs within 7500. 

The observed VDP considering all cluster members {black 
dashed curve) is not well matched by any model profile. 
While the predicted VDP considering all “spectroscopic” 
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Figure 19. Predicted velocity dispersion profiles as a function of accretion 
epoch for model galaxies in our stacked sample of 75 clusters from the Mil¬ 
lennium simulation at z=0.21. The thick blue curve shows the predicted LOS 
velocity dispersion profile of “spectroscopic” cluster members yet to pass 
within r 2 oo and be accreted into the cluster DM halo. The remaining solid 
curves show the impact on the velocity dispersion of progressively includ¬ 
ing those galaxies accreted into the clusters within the last At Gyr as indi¬ 
cated by the color scale. The black solid curve shows the velocity dispersion 
profile obtained by including all cluster members. Blue symbols with er¬ 
ror bars show the observed velocity dispersion profile of the 24/rm-detected 
cluster members from our stacked sample of 30 clusters. The black dashed 
curve shows the corresponding profile considering all Ad>2x IO'^Mq clus¬ 
ter members. 

cluster galaxies peaks at ~0.2r5oo and drops off sharply to¬ 
wards the cluster core {dark red curve), the observed VDP 
shows no corresponding dip in the cluster core. The aios{r) 
instead declines steadily from its peak value LlSa^ in the 
innermost radial bin, falling to values of ~0.75 over the 
range 1.6-2.8 7500, significantly below the velocity dispersions 
predicted by simulations at these radii. Our observed VD P is 
qualitatively similar to that obtained bv iRines et al.l (120031) by 
stacking the member galaxies of eight z<0.05 X-ray luminous 
clusters; their aios{r) also drops from l.lCT,y in the cluster core 
to 0.8(7;^ by 7200, albeit with marginal evidence for a decline 
within 0.1 7200- At larger radii (>27500), their aiosir) drops to 
values of ^0.5 which is even lower than ours and hence 
poses further problems for the simulations. 

One possible explanation is that the mis-match is linked to 
the prediction of too many model galaxies in the cluster core 
{rproj 0.1 7500; Fig- ID, which assuming that the excess pop¬ 
ulation were all accreted early, could artificially increase the 
contribution from low-velocity virialized cluster members, re¬ 
ducing the iJi, estimates for each cluster, and pushing the re¬ 
sultant model curves upwards. 

Irrespective of the difficulties in reproducing the observed 
VDP of all cluster members, the key finding that the veloc¬ 
ity dispersion of star-forming galaxies is 10-35% higher than 


0.0-0.5r„„ 0.5-1.0r,.„ 



Figure 20. Stacked relative LOS velocity distributions of “spectroscopic” 
cluster galaxies for the 75 most massive clusters in the Millennium sim¬ 
ulation at z=0.21 split into three dynamic sub-populations: (i) infalling 
“spectroscopic” cluster members yet to reach pericenter {blue striped his¬ 
togram)', (ii) back-splash galaxies which have passed through pericenter, are 
on their way back out of the cluster i^radiai>^^ are yet to reach apocenter 
{green striped histogram)', and (iii) the virialized cluster population which 
has passed through apocenter {red solid histogram). Each panel represents a 
different bin in cluster-centric radius as in Fig. [7] The standard deviations (in 
units of ajy) and kuitosis values of each distribution are indicated. 

that of the overall cluster population at all radii, along with 
the apparent sharp peak in the VDP at ^0.3r5oo, unambigu¬ 
ously identifies the star-forming cluster galaxy population as 
recent arrivals. Considering a simple kinematical treatment of 
infalling and virialized cluster galaxies in a cluster-scale grav¬ 
itational potential well leads to \T/V\^l for infalling galax¬ 
ies and \T/V\^l/2 for the virialized population, where T and 
V are the kinematic and potential energies. Thus, the veloc¬ 
ity dispersions of the two populations are naively related by 
O'mfall^ V2a 

virial (IColless & DunnlFI996^ . 

From the first dynamical studies of cluster galaxies, the 
velocity dispersions of spiral galaxies hav e been found to 
be systematically highe r than early types (iTammannI 1 1 972I : 
iMoss & DickensI 119771) . Based on much larger samples, 
the stacked velocity dispersions of blue/emission-line galax¬ 
ies were found to be 20% higher t h an the remaining in ¬ 
active galaxies jBiviano et al.l ll 9971: lAguerri et al.l 120071) . 
iBiviano & Katgertl (l2004h showed that if early-type galax¬ 
ies are assumed to have isotropic orbits within clusters, as 
supported by their Gaussian velocity distributions, the kine¬ 
matic properties of late-type spirals are inconsistent with be¬ 
ing isotropic at the > 99% level. Instead they indicate that 
spirals and emission-line galaxies follow radial orbits in clus¬ 
ters, pointing towards many of them being on their first cluster 
infall. 

Figure |7] showed that the LOS velocity distribution of star¬ 
forming cluster galaxies within 7500 to have a rather flat, top- 
hat profile, a high LOS velocity dispersion and a negative 
kurtosis (72"^-0.81) strongly inconsistent with the Gaussian 
distribution typical of a virialized cluster population. This 
flat-topped distribution is well reproduced by model cluster 
galaxy populations that are either infalling into the cluster for 
the first time, or back-splash galaxies which are currently re¬ 
bound^ out of the cluster and are yet to reach apocenter 
(Fig. l20b . The velocity distributions of these two dynami¬ 
cal sub-populations appear indistinguishable within 7500 , both 
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Figure 21. Likelihood function of the quenching time-scale, tg, based on 
fitting the distribution of star-forming galaxies in the caustic diagram, nor¬ 
malized so that f C(tQ) = 1 with tg in units of Gyr. The best-fit value of tg is 
indicated by the vertical dot-dashed line, while the shaded regions indicate the 
1, 2 and 3-f7 confidence limits in tg. The magenta dashed line a nd hashed re- 
gion indicates the best-fit to value and l-cr confidence limits bv IHaines et al.l 
(Ml) based on the systematically low specific-SFRs of stai'-forming cluster 
galaxies within r 2 oo- 


having velocity dispersions ^l.2a^ and negativekurtosis val¬ 
ues 72~-0.8. At l-2r5oo the velocity distribution of star¬ 
forming galaxies becomes more rounded, albeit still with a 
negative kurtosis (72^-0.6), which again is well reproduced 
by the model infalling galaxy population (72=-0.59). At 
these radii, the back-splash population is expected to show 
a more Gaussian-like distribution with 72=-0.09, inconsis¬ 
tent with observations. However, there are only expected 
to be 40% as many back-splash galaxies in this radial bin 
as infalling ones, and so we cannot rule out the possibility 
these star-forming galaxies represent a mixture of infalling 
and back-splash populations. 


5.4. Distribution of galaxies in the caustic diagram 

The information gained from the galaxy surface density 
profiles and VDPs can be combined by comparing the ob¬ 
served spatial distribution of star-forming galaxies in the 
stacked caustic diagram (Fig.|6]l with those obtained from the 
Millennium simulation (Fig.lOll. Using the stacked sample of 
model galaxies from the 75 clusters extracted from the Millen¬ 
nium simulation, the distribution of model star-forming galax¬ 
ies in the phase-space diagram ^vios/ctu versus rproj/rsoo was 
determined using the adaptive kernel estimator, for a range of 
quenching time- scale s from 0-10 Gyr. 

As in Section 15.21 model star-forming cluster galaxies are 
given initial SFRs taken at random from our observed sam¬ 
ple of coeval field star-forming galaxies (SFR>2M0yr“'; 
0.15<z<0.30). These are then set to decline exponentially on 
a quenching time fg once they pass within r2oo of the cluster, 
until their SFRs fall below our nominal limit of 2.OM0yr“'. 
The redshift at which each model galaxy was accreted into 
the cluster is known, from which the number of quenching 
time-scales fg passed between the epoch of accretion and the 
epoch of observation (z=0.21) can be determined, and hence 
its final SFR and whether it would still be classified as star 
forming. Those galaxies yet to be accreted are assumed to still 
be identified as star forming. For each quenching time-scale, 
fg, the spatial distribution of star-forming model galaxies in 
phase-space, pirproj/rsoo, Ai/ios/crO, is determined using the 
adaptive Gaussian kernel method. Each model galaxy i is 
represented by a 2D Gaussian kernel of width 
where p is the geometric mean of the local densities pi of 
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Figure 22. Stacked phase-space diagram, Psvzjav versus rprof rim. for 
model star-forming galaxies ai'ound the 75 most massive clusters in the Mil¬ 
lennium simulation at z=0.2l for the best-ht quenching model, in which star 
formation declines exponentially on a time-scale fg=L73Gyr in galaxies 
from the moment they pass within r 2 oo of the cluster for the first time. Each 
point marks a Af >2x 10*® Mq star-forming galaxy with SFR>2.OM0yr“*, 
colored according to time elapsed since it was accreted into the cluster (mid¬ 
blue if it is yet to be accreted). The grayscale contours indicate the resulting 
probability distribution function of the model stai'-forming galaxies in phase- 
space, P(rproy/r 5 oo, Ai/z/cti/). Each contour indicates a factor 0.2dex change 
in phase-space density of model stai-forming galaxies. The thick contours 
indicate a factor 10 change in phase-space density. 


the model galaxies in phase-space, and ao is the initial ker¬ 
nel width 0.2. To remove the effects of the discontinuity at 
rproj=0, the phase-space distribution is mirrored about both 
velocity and radial axes. The distribution is normalized to 
unity when summed over the region 0.0<(rp,oj/r5oo)<3.2, so 
that it can be considered a probability distribution function, 
Av[qsI of}. 

The probability distributions of model star-forming galax¬ 
ies for each value of fg are then compared with the ob¬ 
served distribution of star-forming galaxies in the stacked 
phase-space diagram (Fig. |6l). The best-fitting model clus¬ 
ter population is identified using a maximum-likelihood anal¬ 
ysis, determining the value of fg for which the likelihood 

A=Y^^^iP{ri,projJrsm,Avilov\tQ) is maximized, taking into 
account corrections for spectroscopic incompleteness and ra¬ 
dial variation in coverage by our 24/im images as before. 

Figure |2T] displays the resulting likelihood function £(fg) 
as well as the 1, 2 and 3 -(t confidence limits in fg. The closest 
match to the observed distribution of star-forming galaxies is 
obtained for a value of fg=1.73±0.25Gyr. Quenching time- 
scales below 1 Gyr and above 3 Gyr are both excluded at >3 cr 
level, primarily due to the radial distribution of star-forming 
galaxies observed in our clusters (Fig.fTSl). 

Figure|22]shows the corresponding stacked caustic diagram 
of model star-forming cluster galaxies for fg=1.73 Gyr. Each 
cluster galaxy is color coded according to its accretion epoch, 
and the overall phase-space density distribution is shown 
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by the grayscale contours. The model star-forming galax¬ 
ies are most frequently found at large cluster-centric radii 
{fproj^i -5 rsoo) and low velocity offsets ( \Ai/[os |1 • 1 cr^). The 
phase-space density of star-forming galaxies drops by a fac¬ 
tor ^3x as rproj/fsoo 0 , for |Ai//oj|<1.4CTi„ producing the 
vertical contours in the left-center of the plot. 

Averaging over the 75 simulated clusters, 42% of the model 
star-forming galaxies with rproj<r 2 oo and line-of-sight ve¬ 
locities identifying them as “spectroscopic” cluster members 
are physically located outside r 2 oo at the time of observa¬ 
tion (z= 0 . 21 ) and have not been inside this radius at any time 
{mid-bluepoints), 33% have been accreted within the last Gyr 
{light blue/green points), 19% were accreted 1-3 Gyr ago, 
and just 6 % were accreted >3 Gyr prior to observation. With 
the model SFRs declining exponentially upon accretion on a 
time-scale of 1.73 Gyr, overall the model SFRs and specific- 
SFRs of the star-forming cluster galaxies {rproj<r 20 o) are re¬ 
duced by 25% from their values prior to accretion, consistent 
with the systematic reductio n of 28% in their specific-SFRs 
found in iHaines et al.l (120131) . 

6. DISCUSSION 

We have shown via two independent methods that star for¬ 
mation in galaxies infalling into clusters is not extinguished 
immediately upon their arrival into the cluster, but requires a 
significant time of the order 1-3 Gyr to be quenched. First 
the steadily increasing surface density of star-forming galax¬ 
ies towards the cluster core (Fig. [TSI l. as opposed to the flat 
radial profile predicted for instantaneous quenching models, 
simply requires star-forming galaxies to survive for a certain 
period within the cluster to build up the over-density seen in 
the cluster core, with a best-fit exponential quenching time- 
scale fg=2.2±0.4Gyr. Second, the velocity dispersion pro¬ 
file of star-forming cluster galaxies is consistently 10-35% 
higher than inactive cluster galaxies at all radii, rising up to 
a sharp peak of 1.44ai, at 0 . 3 r 5 oo, which is inconsistent with 
instantaneous quenching models, but similar to the predicted 
VDPs of models in which star-forming galaxies survive for 
0.5-2.2 Gyr after being accreted (Fig. l20li. Combining both 
radial and velocity information, we compared the spatial dis¬ 
tribution of star-forming galaxies within the caustic diagram 
with predictions from the Millennium simulation, to obtain a 
best-fit value of fg=l .73±0.25 Gyr. 

While these results robustly confirm that star-forming 
galaxies are able to continue forming stars for some signifi¬ 
cant period after being accreted into massive clusters, it is not 
possible simply based upon the kinematics or distribution of 
star-forming galaxies within clusters to distinguish between 
slow quenching models whereby star-formation declines ex¬ 
ponentially (or linearly) over a long time-scale, or a “delayed- 
then-rapid” quenching model in which recent arrivals con¬ 
tinue to form stars normally for a cert ain period At, before 
suddenly stopping (IWetzel et al.ll 2012 h . as demonstrated by 
the similarity of the model curyes in Figures [IT] and [18] 

6.1. The slow quenching of star formation in cluster galaxies 

The key distinguishing feature of the “slow quenching” 
model is its impact on the distribution of SFRs or specific- 
SFRs among star-forming cluster galaxies, systematically 
lowering the mean specific-SFRs as a significant fraction of 
star-forming cluster galaxies are obseryed during this process 
of slow quenching, while the “delayed-then-rapid” quench¬ 
ing model leayes the specific-SFR distribution of star-forming 
cluster galaxies unchanged. Figure [ 8 ] supports the slow 


quenching model by finding kinematic segregation between 
star-forming cluster galaxies with normal or enhanced star 
formation, and those with reduced star-formation, indicatiye 
of ongoing quenching. This suggests that the process of 
quenching occurs oyer a sufficiently long time-scale that the 
kinematics and cluster-centric radii of quenching star-forming 
galaxies to haye ey olyed significantly. More definitiyely, in 
iHaines et al.l (120131) we found that the specific-SFRs of mas- 
siye (^IO'^Mq) star-forming cluster galaxies within r 2 oo to 
be systematically 28% lower than their counterparts in the 
field at fixed stellar mass and redshift, a difference signifi¬ 
cant at the 8 . 7(7 leyel. The entire specific-SFR distribution 
was seen to be shifted to lower yalues, marking the unam¬ 
biguous signature of star formation in most (and possibly all) 
star-forming galaxies being slowly quenched upon their ar- 
riyal into massiye clusters. Assuming a model in which the 
SFRs decline exponentially upon passing within r 2 oo, we ob¬ 
tained a best-fit quenching time-scale of 1.1 "^45 Gyr (ma¬ 
genta dashed line and hashed region in Fig. 12lb . consistent 
with the time-scales obtained here. 

As well as skewing the specific-SFR distribution of cluster 
galaxies on or near the star-forming main sequence, the slow 
quenching model is expected to result in significant numbers 
of galaxies with specific-SFRs of ^10“"yr~* (or SFRs of 0.1- 
1 Mq yr“'), well below the sensitiyity of our Spitzer data, and 
filling in the “ green yalley” gap bet ween star-forming and pas- 
siye galaxies. IWetzel et alJ (120131) found no eyidence of this 
large transition population, leading them to prefer a delayed- 
then-rapid quenching model. One possible way of reconcil¬ 
ing these two results would be a slow-then-rapid quenching 
model, whereby star formation is slo wly quenched for th e 
first 1-3 Gyr, to explain the results of IHaines et alJ (120131) . 
followed-by a second short phase in which the residual star- 
formation is rapidly terminated, in order to retain the obseryed 
bimodal specific-SFR distribution of sat ellite/cluster galaxies 
(IHaines et alJl2()llbl[W^el et al.ll2013h . 

iTaranu et al.l (12014 ) studied the bulge and disk colors of gi¬ 
ant galaxies in z< 0.1 clusters, finding shallow, gradual radial 
trends in disk colors that could be reproduced by slow quench¬ 
ing models similar to our own, but requiring slightly longer 
time-scales of fe=3-3.5Gyr. They ruled out short (<1 Gyr) 
quenching time-scales and “delayed-then-rapid” quenching 
models as both produced much larger and sharper radial 
changes in the median disk colors than obseryed. A number of 
other studies haye also argued for relatiyely long time-scales 
(^l^Gyr) for the quenc hing of star formation in recently 
accreted clu s ter ga l axies (iBalogh. Nayarro & MorrisI 120001: 
Moran et'di l2007t Finn et al. 200i^~ yon der Linden et alJ 


2()10tlDe Lucia et al.ll2012l:IWetzel et al.ll2013l). 

iMuzzin et al.l (120141) found a population of post-starburst 
galaxies in nine z~l clusters, whose stacked spectrum could 
be fit by the rapid quenching (fQ=0.4lg 4 Gyr) of typical star¬ 
forming galaxies, and which traced a coherent “ring” at 0.25- 
0.50 7200 in the stacked caustic diagram that could be repro¬ 
duced by recently accreted galaxies quenched 0.1-0.5 Gyr af¬ 
ter passing within 0.5 7200 - This much more rapid quench¬ 
ing in high redshift clusters could be due t o the shorter gas 
consu mption time-scales of galaxies at (ICarilli & Waited 
120131) or ^ 1.7 X shorter cluster crossing time-scales. 

Slow quenching on ^2 Gyr time-scales matches predic¬ 
tions of staryation models, in which infalling galaxies are 
stripped of their diffuse gaseous halos as they pass through 
the ICM, preyenting further gas accretion onto the galax- 
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ies from the surrounding inter-galactic medium (iLarson et al.l 
[T^ iBekki etaH 12001 iMcCarthv et al.ll2009h . The galaxy 
then slowly uses up its existing molecular and Hi gas reser¬ 
voir over a period of 2-3 Gyr, based on the gas c onsumption 
time-scales observed for nearby spiral galaxies (iBigiel et al.l 
1201 UlBoselli et alll2014l) . This process may be effective well 
beyond the virial radius, with the extended gaseous halo of 
clusters remaining hot (>10®K) and sufficiently dense to strip 
the hot gas atmospheres of infalling g alaxies out to ~ 5 r 2 oo 
(iBahe et allKllBI: iGabor & D^l2()fl . 

Moreover, the growth of galactic DM halos is suppressed by 
tidal effects due to the presence of nearby cluster-mass halos 
(iHahn et al.ll2009h . which results in the peak mass of galac¬ 
tic halos inf alling into clusters oc curring at ^1 .8 r,,,> (3.5 rsoo) 
on average (iBehroozi et al.ll20T^ . Given the tight correspon- 
dance between DM mass accretion and gas accretion onto 
galactic halos, including a significant fraction of gas accreted 
onto t he galaxy via cold, dense filamentary streams at all red- 
shifts (Ivan de Voort et al1l201 Ih . the radius of peak halo mass 
marks the end of continual gas replenishment of the ISM, sig¬ 
nalling the beginning of the end for s tar formation in the host 
galaxy (ISanchez Almeida et al.ll2014h . 

As these galaxies continue their journey towards the clus¬ 
ter core, the ICM they encounter becomes increasingly 
dense and the resultant ram pressures {PramOf-picM^^) become 
strong eno ugh to progressively strip th eir gas disks from the 
outside in (iBriiggen & de Lu^ |2008| ). producing truncated 
Hi, H? gas and Ha disk s (e.g. iKoopmann & Kenne\^l2004l 
iBoselli et al.ll200d 120141) . and outer regions showing re cently 
quenched stellar populations (ICrowl & KennevI 120081) . Hy- 
drodynamical simulations predict that the moderate ram pres¬ 
sures acting on infalling massive spirals as they pass within 
- T 200 are sufficient t o strip half their gas cont ents on 
500-1000Myr time-scales (iRoediger & HenslCTll2005h . while 
observations confirm that ram pressure stripping is effective at 
removing gas and quenching star format ion as far out as rsoo 
(IChung et al.ll20()7tlMerluzzi et al.ll2013l) . The effective time- 
scale for quenching via ram-pressure stripping then becomes 
the ^ 1 Gyr required for galaxies to travel from the cluster out- 
skirts to the pericenter wher e the peak in ram pressure occurs 
(IRoediger & Hensl^ 120051) . The complete absence of Hl- 
nor mal spirals within the rsno radius of Virgo and Coma clus¬ 
ters (IBoselli & Gavazzill200(^ . and the ~0.5-0.8 dex Hi defi¬ 
ciencies of Virgo spirals at fixed stellar mass, NUV - r color 
and stellar mass surface density, cannot be reproduced by 
predictions from chemospectrophotometric models involving 
starvation alone and can only be e xplained by ram-pre ssure 
stripping actively removing the gas (ICortese et al.ll 20 lTl) . 

In typical spirals, the dust-to-gas rati o and internal extinc¬ 
tion (A fuv) decline steadily with radius (iMunoz-Mateos et al.l 
Hool), meaning that for ram pressure stripping events where 
gas is removed from the outside in, the more extended un¬ 
obscured star formation component should be preferentially 
quenched prior to the more concentrated and bound obscured 
SF. This is supported by the finding of a significant cluster 
population of 24/rm-detected spirals with reduced SFRs, in¬ 
dicative of ongoing slow quenchingj redd ened by dust onto 
the optical red sequence (IWolf et al.ll2009h . This appears in¬ 
consistent with our finding that the quenching time-scale for 
UV-selected star-forming cluster galaxies is ^IGyr longer 
than that for 24/rm-selected galaxies. We suggest that the 
longer Iq for the NUV-selected star-forming galaxies is due 
to the 1-1.5 Gyr time-scale required for recently quenched 


galaxies to migrate from the blue clou d to red sequence i n 
the (NUV-r) versus C-M diagram (iKavirai et al.l 120071) . 
while the 24/im emission mostly originates from Hll regions 
with ongoing star formation (e.g. lCalzetti et al.|l2005l) . 

6.2. Pre-processing 

In the SF-radius trends of Figs. |2]-0] the fraction of star¬ 
forming cluster galaxies rises steadily with radius, but cru¬ 
cially the fspif) remain stubbornly 20-30% below that seen 
in coeval field populations even out at 5 rsoo, well beyond the 
maximal distances back-splash galaxies are expected to re¬ 
bound to. As a result, these trends cannot be reproduced by 
models in which star formation is quenched in infalling galax¬ 
ies via processes that are only initiated when they are accreted 
into the clusters (Fig. [Tsll. Our SF-r adius trends are very sim¬ 
ilar to those of IChung et al.l (1201 ll) . who analysed 69 z<0.1 
clusters covered by both SDSS spectroscopy and WISE 22/rm 
photometry, finding that the fraction of star-forming galaxies 
with L/«>4.7x 10 '°Lq increases steadily with cluster-centric 
radius, but remains well belo w the field value even at ~3 r 2 oo- 
Ivon der Linden et'dl (1201 011 found that suppression of star- 
formation in cluster galaxies could be traced out to ~ 4 r 2 oo- 
IWetzel et al.l (1201 2^ found their fspir) suppressed with re¬ 
spect to field values as far out as 10 r 2 oo around cluster-mass 
halos (M 2 oo> 1 O*"^M 0 ), but interestingly saw this suppression 
entirely limited to satellite galaxies, with no evidence of sup¬ 
pression among centrals found beyond 1-2 r 2 oo- 

The best way to reproduce the observed SF-radius trends 
appears to require a certain fraction of infalling galaxies 
to arrive onto the clusters having already been quenched 
(Fig. M- One commonly identified mechanism by which 
galaxies may be transformed at large distances from the clus¬ 
ter center is through “pre-processing” in infalling galaxy 


groups (e.g. iKodamaet al.lF2 001l: iFu iital 12004 iBerrier et~ffl 
l2009i iMcGee et al.l 120091: iDressler et al.l 12013 ). Star forma¬ 


tion is suppressed in group galaxies, with fsp values inter- 


tions (iWilman et al.l2005HBai et al.ll2010l; 

McGee et al.||201 1 

Rasmussen et al.l2012^ up to at least z~l ( 

Balogh et al.KOll 

Ziparo et all 120141). The fraction of star-forming galaxies 


within groups declines steadily with increasing group mass 
(at fixed stellar mass an d group-centric distance) and prox¬ 
imity to the group cente r (IWeinmann et al.ll2006l: IWetzel et al.l 
1201:^ IWoo et aITl2013l) . The fractions and specific-SFRs of 
star-forming galaxies within groups both show accelerated 
declines since 1 w ith respect to the coeval field population 
(iPopesso et al.ll2()T4h and also “filament-like” environments 
with comparably high g alaxy densities but no X-ray emis¬ 
sion (IZiparo et al.ll2014ft. indicatiye of o ngoing slow quench¬ 
ing within groups (fBalogh et'dl 1201 ll). comparable to that 
seen in our clusters iHaines et al.ll2013l) . Galaxies in groups 
with masses > 1 O'^M 0 are Hl-deficient at fixed stellar mass 
and NUV - r color, suggesting that ram-pressure stripping 
can remoye atomic gas from spirals in such groups 
(ICatinella et al.ll2013h . 

Galaxy groups are ubiquitous, and hos t ^50% of galaxies 
in the local Uniyerse (iMcGee et al.ll2009h . The impact of en- 
yironmental quenching in groups will be manifest also in a 
reduced global fraction of star-forming galaxies in our coeyal 
field sample. Hence, in order to reproduce the obseryation that 
the fraction of star-forming galaxies in the infall regions of 
clusters (2-3 r 2 oo) remains significantly lower than that seen 
in coeyal field galaxies, a specific form of “pre-processing” 
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Figure 23. History bias for galaxies falling into massive clusters. The solid 
curves show the fraction of model galaxies surrounding massive clusters 
(Mj./> 1 O*^'^M 0 ) at z= 0.21 which already reside in halos over a given mass 
as a function of projected cluster-centric radius. Only those galaxies yet to 
pass within r 2 oo and be accreted into the cluster are included. The dashed 
curves show the corresponding cosmic fractions for field galaxies over the 
whole Millennium simulation. 

is required, in which the galaxies infalling into clusters are 
more likely to be existing members of groups and hence more 
likely to be “pre-processed” than the cosmic average. That is 
we require the mass function of DM halos hosting galaxies in 
the surroundings of clusters (not including the cluster itself) 
to be top heavy and biased towards group-scale masses with 
respect to the cosmic av erage, as expected in suc h over-dense 
regions of the universe (iFaltenbacher et al1l2010h . 

To examine this issue, we construct the mass function of 
the DM halos hosting galaxies at different cluster-centric radii 
and compare it with that of coeval field galaxies. We se¬ 
lect galaxies within 4000 km s“* of the mean redshifts of rich 
clusters ( Mhain >1O*^'^M0 ) in the Millennium simulation at 
z=0.21 (see lLu et al.l2012l for details), excluding those galax¬ 
ies already accreted into the clusters. Figure 123 shows the re¬ 
sulting fractions of these model galaxies which reside in halos 
over a given mass as a function of projected cluster-centric 
radius (solid curves) for five values of minimum halo mass 
from lO'^M© to lO'^^M©. The dashed lines show the cor¬ 
responding fractions for coeval field galaxies over the whole 
Millennium simulation. This confirms that galaxies infalling 
into clusters are indeed more likely to already reside within 
DM halos of a given mass than coeval field galaxies, for each 
halo mass range. While the effect is marginal for lO'^M© 
mass halos, galaxies infalling into massive clusters are ~35% 
more likely to already be residing in group-mass halos with 
A/haio>10^^MQ and twice as likely to already be residing in 
4^haio>10*"^M© halos than their counterparts in the field. 

We find numerous potential sites where infalling galaxies 
can be pre-p rocessed before bei ng accreted into the LoCuSS 
clusters (e.g. lPereira et al.lbOlOl) . Across the 23 clusters with 
XMM imaging, a total of 30 X-ray groups (with extended X- 
ray emission detected at >4(7 levels) have been identified with 
redshifts placing them inside the cluster caustics, indicating 
that they are most likely infalling into the primary cluster 
(Haines et al. 2015, in preparation). Only six further “iso¬ 
lated” X-ray groups were detected (at >4 (t) in the same XMM 


images over the rest of the 0.15-0.30 redshift range. 

The observed ratio of ~5 X-ray groups associated with the 
clusters for every field X-ray group in the remainder of the 
0.15<z<0.30 volume covered by XMM, is double (quadruple) 
the ratio obtained for cluster and field Mk<M’^+\.5 (24/im- 
detected) galaxies in the same two volumes. This suggests 
that infalling galaxies are 2-Ax more likely to be members of 
X-ray groups than in typical field regions, assuming that both 
group samples have similar numbers of members per group. 

A second factor in the shortfall in fsp at large radii could 
arise from a bias in the mass assembly history of the in¬ 
falling galaxies and groups themselves. iMaulbetsch et al.l 
(l2007h find that galactic ~10'^M© halos in high-density re¬ 
gions, such as the infall regions of massive clusters, form ear¬ 
lier, have more active merger histories, and have much lower 
mass accretion rates (an d commensurate gas accretion rates; 
Ivan de Voort et al.l 1201 Ih at late epochs (z<0.5) than those 
which form in low-density field environments. Such halos 
in high-density regions are ^4x likelier to not be accreting 
any mass (or even losing it) at late epochs than those in low- 
density regions. 

7. SUMMARY 

We present an analysis of the radial distribution and 
kinematics of star-forming galaxies in 30 massive clusters 
at 0.15<z<0.30, combining wide-held Spitzer 24pm and 
GALEX NUV photometry with highly-complete spectroscopy 
of cluster members. To gain insights into how the observed 
trends relate to the continual accretion of star-forming spirals 
onto massive clusters and subsequent quenching of star for¬ 
mation, we follow the infall and orbits of galaxies in the vicin¬ 
ity of the 75 most massive clusters in the Millennium cos¬ 
mological simulation, obtaining a series of predicted model 
trends that should have general applicability for understand¬ 
ing galaxy evolution in cluster environments. Our main re¬ 
sults are summarized below: 

1. The surface density of star-forming galaxies declines 
steadily with radius, falling ^ 15 x from the cluster core 
to 2 7200- This simple observation requires star for¬ 
mation to survive within recently accreted spirals for 
2-3 Gyr to build up the apparent over-density of star¬ 
forming galaxies within clusters. 

2. The velocity dispersion prohle of the star-forming clus¬ 
ter galaxy population shows a sharp peak of 1.44 at 
0.3 7500, and is 10-35% higher than that of the inac¬ 
tive cluster members at all cluster-centric radii, while 
their velocity distribution shows a flat, top-hat profile 
within 7500- All of these results are consistent with star¬ 
forming cluster galaxies being an infalling population, 
but one that must also survive ~0.5-2 Gyr beyond pass¬ 
ing within 7200 to achieve the high observed velocities. 

3. The distribution of star-forming galaxies in the stacked 
caustic diagram are best-fit by models in which their 
SFRs decline exponentially on quenching time-scales 
fQ=1.73±0.25 Gyr upon accretion into the cluster. The 
above results, and the observed kinematic segregation 
of star-forming galaxies accordi ng to their specific- 
SFRs, support the conclusion from lHaines et al.l (120131) 
that star formation in most (and possibly all) high-mass 
star-forming galaxies is slowly quenched on accretion 
into rich clusters on 0.7-2.0Gyr time-scales 
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4. The fraction (fsp) of star-forming cluster galaxies rises 
steadily with cluster-centric radius, increasing five-fold 
by 2 r 2 oo, but remains well below field values even at 
3 r2oo- Pre-processing in infalling galaxy groups ap¬ 
pears the most likely explanation for this suppression 
of star-formation at large distances from the cluster. 
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